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We present a self-contained theoretical and computational framework for dynamics following 
photoexcitation in quantum dots near planar interfaces. A microscopic Hamiltonian parameterized 
by first principles calculations is merged with a reduced density matrix formalism that allows for 
the prediction of time-dependent charge and energy transfer processes between the quantum dot 
and the electrode. While treating charge and energy transfer processes on an equal footing, the 
non-perturbative effects of sudden charge transitions on the Fermi sea of the electrode are included. 
We illustrate the formalism with calculations of an InAs quantum dot coupled to the Shockley state 
on an Au[lll] surface, and use it to concretely discuss the wide range of kinetics possible in these 
systems and their implications for photovoltaic systems and tunnel junction devices. We discuss the 
utility of this framework for the analysis of recent experiments. 



I. INTRODUCTION 

Nanostructured materials represent one of the most 
promising routes for the creation of novel energy harvest- 
ing and optoelectronic devices. One of the key challenges 
in this field is navigating the vast design space to search 
for an appropriate combination of material properties for 
a specific application. At present, the solution to this 
fundamental problem can at best be constructed for spe- 
cific classes of systems. In the present paper, we tackle 
this problem for one of the key configurations that has 
emerged in the fields of photovoltaics and nanoelectron- 
ics: quantum dots acting as chromophores in the vicinity 
of semiconductor and metallic surfaces 1-3 . 

We focus on charge kinetics in these systems and con- 
struct a microscopic dynamical theory that is capable 
of both describing the observable experimental phenom- 
ena, and of making quantitative predictions for future 
research. The key challenge addressed in the paper is 
the construction of a framework in which a strong con- 
nection is maintained between the microscopic param- 
eterization of the Hamiltonian, and the description of 
dynamics within a restricted set of excited states. Thus 
key insights from an efficient navigation of the parameter 
space can be related to specific physical properties of the 
underlying materials in these systems. 

Theories based on model Hamiltonians 4-8 allow one to 
construct a picture of the possible regimes of dynamics, 
but the connection to realistic modeling of specific ma- 
terials may be hard to establish. On the other hand, 
parameter free ab-initio studies 9,10 can only address rel- 
atively small systems and the key physical features of 
quantum dot systems can not be treated. The work pre- 
sented in this paper may be viewed as a compromise be- 
tween these two extremes. Starting from the fundamental 



many-body Hamiltonian for the light-matter interaction, 
we derive a low energy effective Hamiltonian in which 
all matrix elements are calculated from the single parti- 
cle wavefunctions of the subsystems, each represented by 
a well-established model. The wavefunctions and other 
characteristics of the surface are obtained from an ab- 
initio Density Functional Theory based calculation. The 
frontier states in the semiconductor quantum dot are ob- 
tained from established effective mass models. In the 
region of overlap between these subsystems, these wave- 
functions can be treated on equal footing. In this way, 
our low energy effective Hamiltonian is fully derived, with 
the addition of no parameters. This Hamiltonian is then 
used to develop a master equation for the reduced dynam- 
ics of the quantum dot. Simulations based on the mas- 
ter equation provide quantitative analysis of the range 
of possible charge kinetics in these systems with a clear 
connection to the constituent surface and quantum dot 
materials. 

By constructing the Hamiltonian in this way, a mean- 
ingful analysis of charge and energy transfer channels in 
isolation and in competition with each other emerges. 
Thus one can pursue important questions about charge 
kinetics such as: How does the extinction of optical power 
by the quantum dot change in the presence of a surface? 
How is cooling of hot carriers affected by the presence 
of energy transfer to the surface? How does electro- 
static coupling of the quantum dot to a surface affect the 
lifetime of exciton, bi-exciton, and other multi-electron 
states? How does non-radiative recombination of exci- 
tons in these systems affect current extraction and pho- 
tolumincsccnce? Our framework is capable of answering 
these questions concretely. We illustrate this within a 
specific system of InAs quantum dot on a gold surface. 

In recent years, several studies of kinetic processes in 
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models for single molecule junctions under the influence 
of both applied electrical bias and optical fields have been 
published 4-13 . The studies of asymmetric dipole coupling 
in the steady state conductance 5 ' 8 have revealed photo- 
induced current generation, and current induced photo- 
emission. This scenario, however, is inapplicable to single 
interface systems in which only a photo-generated cur- 
rent can can exist, and interfacial polarization plays a 
different role, as described below. Electron dynamics 
in molecular chromophores at semiconductor interfaces 
has been studied in depth for a small number of systems 
and ideal situations fully ab-initio 10 . While these studies 
are invaluable for quantitatively settling many questions 
about the microscopic processes and their dynamical in- 
terplay, they do not capture the aspects arising from the 
relatively large sizes of quantum dots. Furthermore, the 
treatment of image potential at the interface remains an 
external input to any computation based purely on den- 
sity functional theory 14-16 . 

The highly polarizable surfaces of the planar elec- 
trode and the quantum dot lead to electrostatic interac- 
tions that significantly affect the quasiparticle and opti- 
cal bandgaps, tunneling rates, as well as energy transfer. 
For a spherical quantum dot, an important fundamental 
effect of the presence of a planar surface is the formation 
of a dipole moment and strong corrections to multipole 
moments of its exciton states. Thus, the polarizability of 
electrode surface becomes a mechanism for non-radiative 
recombination of excitons via energy transfer to the elec- 
trode. In our numerical results, we show the effects of 
this polarizability and quantify the significance of high 
order multipole moments of charge distribution on the 
non-radiative exciton recombination in the narrow gap 
In As quantum dot. Having a microscopic Hamiltonian 
in hand further allows us to compare this exciton decay 
with the dissociation across the junction. This type of 
analysis, for example, is fundamental to optimization of 
current extraction in photovoltaic applications of these 
systems. 

Furthermore, the size of a quantum dot also yields 
many closely spaced energy levels, which lead to dynam- 
ical effects that often do not arise in single molecules, 
especially if a simplified treatment is limited to just the 
highest occupied and lowest unoccupied molecular or- 
bital. The spacing and the number of energy levels qual- 
itatively affects the dynamics of charge injection, energy 
exchange, and the recoil of a hole (electron) in response to 
the tunneling of an electron (hole) to the electrode. Un- 
derstanding the time for build-up of these transitions in 
relation to the magnitude of these transition rates is fun- 
damental to determining the regimes where a Markovian 
description of charge and energy transfer breaks down. 
In experimental terms, it allows us to understand when 
to expect deviations from Lorentzian lineshapes in both 
linear absorption and non-linear optical spectroscopy of 
these systems. The non-Mar kovian effects originate phys- 
ically from the coupling of quantum dot states to the 
electrode, as has also been discussed by Fainberg et al. 4 



in the case of molecular junctions. However, it can be 
affected significantly by the spacing of energy levels in 
the quantum dot, which is a complementary aspect that 
exists naturally in our work. Furthermore, the levels also 
couple significantly, via the Coulomb interaction, to the 
incoherent particle-hole excitations in the Fermi sea of 
the electrode, which opens additional channels of energy 
transfer beyond those discussed in studies of molecular 
systems. 

With a few notable exceptions 4,17 ' 18 , studies in molec- 
ular transport generally treat the electrodes as passive 
Fermionic reservoirs, thus ignoring the scattering of elec- 
trons in the leads as a result of the excitation dependent 
Coulomb potential of the molecule. Under certain con- 
ditions, discussed in this paper, this potential can sig- 
nificantly alter the transport and relaxation in tunneling 
junctions 19-22 , as well as the optical absorption 23 . For 
example, tunneling of an electron out of the quantum 
dot yields sudden transition of its charge, and can cause 
significant dynamical fluctuations in the surface charge 
density in the electrode. Studies of the effects of this kind 
have a long history in X-ray emission and absorption in 
bulk metals, and optical absorption in doped quantum 
wells 24-26 They are well-known to be composed of two 
competing contributions: the Mahan exciton (ME) 27-30 
arising from the attraction of the electron in the Fermi 
sea to the hole, and the Anderson orthogonality catas- 
trophe (AOC) 30,31 arising from the vanishing overlap be- 
tween the initial and the final many-body state. Together 
they define the phenomenon of the Fermi edge singularity 
(FES). A highly relevant example to the systems consid- 
ered in this paper is the recent observation of ME in 
InGaAs/GaAs quantum dot heterostructures 3 . 

The theory formulated in this paper fully accounts for 
FES phenomenon self-consistently alongside charge and 
energy transfer processes. A novel aspect of the geometry 
considered in this paper is that it is ideal for exploring 
FES within the hole bands of a p-doped electrode. This 
has remained unexplored in FES studies in bulk metals 
and quantum wells because the core charge in this case is 
the much lighter electron, the motion of which diminishes 
the FES signature. On the other hand, the electron lo- 
calized inside the dot presents no such problem and new 
effects arising from scattering in non-parabolic bands and 
the much larger sub-band mixing in hole states can be 
explored. In addition, identifying systems in which these 
effects yield important signatures in optical response and 
tunneling current is important for the correct interpreta- 
tion of experimental data as well as device engineering. 

The FES has also been studied in resonant tunneling 
devices 19 , and was first predicted in these systems by 
Matveev and Larkin 19 . Abanin et al. emphasized the 
tunability of the FES effect by engineering the geometric 
aspects of the system, and elucidated the novel effects 
arising within a non-equilibrium electron gas 21 . The op- 
tical response we formulate here naturally leads to an 
extension of this idea within the interfacial quantum dot 
systems where tuning the relative effect of the ME against 
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the AOC can be achieved by controlling whether the final 
state of the optical excitation lies in the Fermi sea or in 
the capping layer of the heterostructure. 

In a study by Despoja et al 23, on the ab-initio cal- 
culations of the core-hole spectrum of jellium surfaces 
the FES was found to be very weak for a core-hole re- 
siding outside the surface. This is due to very strong 
screening by the free electron gas in the metal. While 
this may be expected for the small overlap between the 
sharply screened potential inside the metal, and the ex- 
tended states of the slab, our numerical calculations show 
it to be true even for the sea of Shockley surface states 
on Au[lll]. On the other hand, we expect the effect 
to increase dramatically for thin films supported on an 
insulating substrate with low dielectric constant. 

Our paper is organized as follows: in Sec. II we derive 
the microscopic Hamiltonian. The detailed expressions 
for all the required Hamiltonian matrix elements are pro- 
vided in Appendix A. In Sec. Ill, we develop our model 
for the dynamics, specifically a time-convolutionless mas- 
ter equation for the density matrix of the quantum dot. 
From this equation, we obtain the dynamical rates of 
charge and energy transfer, as well as the optical response 
including the effects of the FES. In Sec. IV we apply our 
model to calculations of charge and energy transfer for 
an InAs quantum dot on Au[lll] surface, and study the 
effect of the FES as well as the competition between tun- 
neling, cooling, and non-radiative exciton decay on the 
dynamics. In Sec. V we discuss the application of our 
theory to the modeling and analysis of experiments, and 
the possible extension of the theory developed here to 
that of a quantum dot array supporting inter-dot charge 
and energy transport. We also discuss in this section how 
the vibrational modes of the quantum dot, neglected in 
the explicit analysis, can be described within the frame- 
work presented. In Sec. VI we conclude. Details not 
found in the text may be found in Appendices A and B. 



II. MICROSCOPIC HAMILTONIAN 
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We start with the Hamiltonian, 
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where ^(r) is the fcrmionic annihilation field operator. 
The interaction of electrons with light occurs through 
the vector potential A(r,i) associated with the optical 
field. We treat the field as a classical external force in 
this work. Therefore we do not include energy stored in 
this field in the above Hamiltonian. The potentials as- 
sociated with the QD and the electrode are given by the 
functions t/z)(r) and Ul(t) respectively (see also Fig. 1). 
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Figure 1: (a) The geometry of the system studied in 

this paper showing the important parameters for 
determining the electrostatic potentials, (b) Schematic 
profile of potentials and wavefunctions of the electrode 
and the quantum dot in a plane perpendicular to the 
electrode surface and passing through the center of the 

dot. 



In the last term of (1) V(r — r') is the Coulomb interac- 
tion among all fcrmions. We have neglected the phonons 
and electron-phonon interaction in the present model for 
brevity, but it can be incorporated straightforwardly in 
our formalism, and we comment in Section IV on how 
this can be accomplished. 

Our approach is akin to the Bardeen approach to tun- 
neling 32 . We make a physically reasonable distinction 
between states of the QD and the electrode, and then 
develop a theory for the charge and energy exchange be- 
tween the two sets of states. It is then convenient to first 
write the field operator as a sum of field operators that 
create/annihilate particles confined to the QD and the 
electrode, 



#(r) = *£>(r) + * L (r). 



(2) 



With the choice of the single particle basis below, how 
the states within the two subsystems can be identified 
will be detailed. It is implicit in the formalism of the 
reduced dynamics of the QD (see Sec. Ill) that the elec- 
trode states, which are traced over, are orthogonal to 
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the QD states included in the dynamics. Any distinction 
based on exact vanishing of the states within appropri- 
ately defined volumes of the two subsystems would not 
yield an orthogonal basis in general 33 . However, with the 
barrier potential equal to several electron volts, the low 
energy states of the QD and those near the Fermi level of 
the electrode, even when calculated in isolation of each 



other, decay exponentially within the barrier with a char- 
acteristic length of about 1-2 A. The QD and electrode 
states most relevant to the dynamics are then orthogonal 
to a good approximation. We exploit this property in the 
calculations, but note that the approximation lies in the 
choice of basis and not in the formalism. 

We substitute (2) into the Hamiltonian (1) above, and 
group the terms as follows, 



Hbl 



dr$! ] L {r) 
dr& D (r) 



1 A V4 , 
2m I i 



A(r,t)\ +U D (r) + U L (r) + j dr'V{r - r') {*^(r')* D (r') + p img (r')} 
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V + eA{r,t)j +U D {r) + U L (r) + J dr'V(r - r')^ L (r')^ L {r') 



2m \ i 



V + eA(r,t) +U D (r) + U L (r) 



Vdl = j dr j dr'^ D {r)y D (r)V(r-r') \^{{r')^ L (r') - p lmg (r')] 



*u(r), (3) 

(4) 

(5) 
(6) 



The term Hp given by (3) describes the QD states 
in the presence of the electrode potential, the Coulomb 
interaction among carriers, and the classical electrostatic 
interaction with the electrode in the form of an image 
potential. The latter is included via the induced density 
(a real- valued function) pi m g{f), which is later subtracted 
in (6). We call the basis states of Ho for A = the 
QD states, and the matrix elements of —im~ 1 hA(r) ■ V 
between these states defines the optical excitation of the 
QD in the presence of an electrode. 

Similarly, in (4) defines electrode states and their 
optical interaction in the presence of the lattice potential 
Ul and the Coulomb interaction among carriers. The sin- 
gle particle states of Hl are calculated with full atomic 
scale detail, after setting A = 0. Then the Coulomb in- 
teraction appearing in Eq. (4) is incorporated into the 
treatment of the states at the level captured by Den- 
sity Functional Theory. In the particular example we 
use for illustration, only the surface bound states play 
a direct role in calculations, while the remaining states 
are kept formally for completeness. In the final reduced 
dynamics of the QD, the sum of the interaction matrix 
elements over these states yields well-defined single and 
two-particle response functions. The actual problem of 
treating the electrode is then reduced to the calculation 
of these response functions in the presence of a surface. 
The surface bound states are calculated by constructing 
the surface Green functions from the Kohn-Sham density 
functional theory, as discussed at length in Section IV. 

Next, the term Hdl given by (5) describes hybridiza- 
tion, including the optically driven excitations in which 
an electron or a hole is excited directly into a final state 
in the electrode. This term thus describes charge trans- 
fer. The last term, Vdl, given by (6), describes energy 



transfer mediated by the Coulomb interaction between 
the QD and the electrode. The advantage of adding 
and subtracting pi mg (r) in the above expressions is two 
fold: the QD states include the image potential non- 
perturbatively, and the dynamical interaction of these 
states with the electrode reduces to fluctuations around 
this density. As we will see in Section IIIC, this natu- 
rally leads to describing the energy transfer in terms of 
a dynamical longitudinal susceptibility of the electrode. 

In the above expansion we neglect the term 
^\{v)U ^r)^ i,{r) that appears in H^, which would re- 
sult in additional renormalization of the electrode states 
in response to the presence of the neutral QD. We disre- 
gard this term because the exponential decay of the elec- 
trode wavefunctions in the barrier, combined with the 
extended nature of these wavefunctions inside the elec- 
trode, makes the effect of this perturbation very small. 
On the other hand, the changes induced by the presence 
of the electrode has a significant effect on the boundary 
conditions for the states localized to the QD, and thus 
the analogous term \l/j : ,(r)/7i(r) v I' _o(r) is included in the 
Hjj term. We remark that the neglected term may be in- 
cluded by calculating the scattering states starting from 
the eigenstates of Hl 29 - This would perhaps be necessary 
for some mesoscopic electrodes lying at a very small dis- 
tance from the QD, or in cases with small barrier heights 
such that the electrode wavefunctions significantly probe 
U D (r). 

We have also neglected the particular exchange inter- 
action where one carrier lies in the QD state, and the 
other in the electrode state. This interaction should also 
be exponentially suppressed for the low energy and well 
confined states of relevance to this work. Related to this 
exchange is the Coulomb-driven tunneling process which 
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is also neglected in comparison to the contribution by sin- 
gle particle kinetic and potential energies as well as opti- 
cal interactions, Hdl- The exchange interaction among 
carriers within each subsystem is implicit in the above 
expressions. 

Let us now turn to specification of the states compris- 
ing jj and an d begin with a set of single particle 
basis states. Let \g) be the ground state of the QD, and 
fn{f) = (r\ n) be the single particle excited states sat- 
isfying, 



2toq 



V 2 + U Deff (r) + Z(r)-e r , 



<Pn(r) = 0. (7) 



Here Uoeff represents the pseudopotential for a single 
electron or a hole added to the neutral ground state of 
the QD, and we have introduced the electrostatic self 
energy of a point charge in the vicinity of polarizable 
surfaces such as that of the QD and the electrode. The 
electrode contribution is given by 34 



E L (r) = J dr'V(r 



r')pi m g(r'). 



(8) 



This is a well-known formulation of an exact one-body 
potential representing the electrostatic energy stored in 
the polarization reaction field that is induced by a unit 
charge on dielectric surfaces 35-37 . 

We label the solutions to the eigenvalue equation (7) 
as electron states, ip e (r), for the addition of a single elec- 
tron to the QD state e above the quasi-particle energy 
gap, and hole states, yj£(r), for the removal of an electron 
from a valence state, h. While the solutions at higher en- 
ergies significantly violate the boundary conditions at the 
electrode surface, they are nonetheless useful in forming 
a convenient basis set for expanding the multi-particle 
states. The present paper discusses only the exciton as 
the multi-particle excitation, and we express its wave- 
function as, 



$ x (r e ,r h ) 



eh 



ip e (r e )(pl(r h )^ eh . a 



(9) 



where the coefficients <5> e h;x are determined from a varia- 
tional calculation including the Coulomb interaction be- 
tween the electron, hole, and their induced surface polar- 
izations (see (3) and Sec. IV). This approach has been 
employed widely in semiconductor optics 38 ' 39 . 

The electrode states diagonalize Hl in (4) for A = 0, 
and we write them as ij) v k = u u k(r)e lk ' r . While it is not 
essential for the theory presented, we have parameter- 
ized these states by a two-dimensional quasi-momentum 
fc, thus assuming that the electrode has a planar surface. 
In a semi-infinite electrode, v may also take continuous 
values within regions of the bulk excitations in the pro- 
jected density of states. In terms of the (non-orthogonal) 
basis set above, (2) can now be written more precisely as, 



tf(r) 



fn(r)c n + ^2 u uk {r)e lK ' r c 

uk 



ikr , 



vk ; 



(10) 



where the first term on the right hand side corresponds to 
^ d and the second to ^ l- The operators c n and c„fc an- 
nihilate particles in states tp n and ip u k respectively. The 
sum over n is truncated to states that are well localized 
on the QD in the sense that their weight in the electrode 
region is negligible. The remaining states, which are high 
in energy and have a significant fraction of their wave- 
function inside the electrode, are then all viewed formally 
as part of the states to be traced over in the reduced dy- 
namics. 

Separating the optical interaction from the A = form 
of the terms in (3)-(6), we now write the Hamiltonian as 



H = H D + H L + H T + V C + + H { [ ] + 



(11) 



in which the first two terms are the Hamiltonians for the 
quantum dot and the electrode states in the absence of 
radiation, 



\ri) e n (n\ , 



uk 



Hl — } , ^vk^ vh c u k- 



(12) 
(13) 



Here the e n are the energies of the QD states, including 
the ground state of the neutral QD, single electron and 
hole states, and the neutral exciton states. The e u u are 
the dispersion of electrode energies with band index v. 
The third term in (11), which arises from (5) represents 
the hybridization between the QD and the electrode, and 
we write it as, 

H T = \e) (g\T e9 + \h) (g\T h9 

+ \h)(x\T hx + \e)(x\T ex + c.c. (14) 

Here we have introduced the electron transfer operators. 
T ab , corresponding to the tunneling induced change in 
the charge state of the quantum dot. These formal oper- 
ators are introduced to simplify the dynamical model in 
Sec. Ill below, and are defined as, 



vk' 



q-eg _ T eg r 

1 — v 

Thx \ ^ rphx t 

uk 

'-rex \ rpex 



(15) 
(16) 

(17) 
(18) 



vk 



The matrix elements T®£ 



represent tunneling ampli- 
tudes, and can be computed from the electronic states 
ip n and tp v k as described in detail in Appendix A. Recall 
that tp h represents a single electron orbital in the valence 
band 40 , such that yf h is the state of the corresponding 
hole. Thus the matrix element T? describes the trans- 
fer of an electron from state h in the valence band of the 
QD to the electrode. We have now defined all the elec- 
tron transfer operators needed for the formalism below 
to describe addition or removal of electrons from the QD. 
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Returning to (11), the fourth term represents the 
Coulomb interaction in (6), which we expand in QD basis 
states, 

V C = I") H Vnm, (19) 
nm 

where the operator V nm acts on the electrode states, and 
is defined as 

Kin' ~ ^ ^ ^ ^ Vn,vk\n' tV'k'CykCy'k' • (20) 
v,v' k,k' 

The matrix elements in this expression follow directly 
from (6), and are defined in full form in (A 14). For 
n 7^ to, the V nm represent quantum fluctuations around 
the classically induced density of the electron gas and 
lead to energy and polarization transfer between the QD 
and the electrode. On the other hand, the diagonal 
terms, V nni cause random fluctuations in the energy level 
of the QD. This may be thought of as a back-action from 
the excitations in the electrode induced by the QD poten- 
tial. As was discussed in the Introduction, this coupling 
results in the FES and AOC phenomena. 

At this stage we let V nm represent all electronic exci- 
tations of the system so that the matrix elements may be 
taken to be the bare Coulomb interaction. In actual cal- 
culations, however, it is more convenient to identify a set 
of elementary excitations of the electrode strongly cou- 
pled to the QD , and then renormalize this coupling by the 
interactions among these excitations, and their interac- 
tions with the weakly coupled excitations. We will see in 
Sec. Ill that this can be achieved essentially by defining a 
frequency dependent dielectric function for the electrode 
surface, in which the relevant interactions are included 
by construction. 

We now turn to the last three terms of (11), which 
describe the interaction of the entire system with the 
external electromagnetic (EM) field. We first write the 
matrix elements of the velocity operator in the standard 
form as, 

Vrm'W = ^A(t)S nn , (21) 

+ i / *KMV^M-^WVC(r)], 

where <f) n may be set to a QD state Lp n or an electrode 
state "i/v The fundamental optical transition is the exci- 
ton, and we define its matrix element as, 

V xg (t) =J2veh(t)$eh;x, (22) 

e,h 

and write the interaction between the EM field and the 
QD as, 

H { D ] (t) = e]T A(t) • v xg (t) \x) (g\ + c.c. (23) 

X 

Similarly the interaction between the electrodes and the 
EM field is given by 

= e^^A(t)-v„ fcyfe -(t)4c i/ , fe -+c.c.(24) 

u,v' k.k' 



An additional light-matter coupling that we do not 

(r) 

consider in much detail here, Hj, , describes the radia- 
tion driven charge transfer between the QD and the elec- 
trode. This term introduces an externally controllable 
exciton dissociation between the QD and the electrode, 
and takes the form, 

H { t \t) = \e) (g\ ^ A(t) • v e;i , fe (t) Ci , fe (25) 

e vk 

k + c - c - 

h i>k 

This completes our construction of the microscopic 
Hamiltonian, and we now turn to the description of the 
charge dynamics of an electrode coupled quantum dot. 

III. DYNAMICS FOLLOWING 
PHOTO-EXCITATION 

We study dynamics within the restricted Hilbert space 
of four classes of states: the ground state \g), single elec- 
tron states, |e), single hole states \h) , and the exciton 
states \x). This restriction is only for convenience, and 
may be lifted by expanding the set to include bi-exciton 
states and even larger multi-electron complexes. Within 
each class, however, we do allow for an arbitrary num- 
ber of states to exist. Since this system is coupled to 
the electrodes, the unitary evolution governed by the 
Schrodinger equation applies to the full density matrix, 
which we denote as £(i), and it obeys the equation, 

|e(*) = ^[ff (*),€(*)], (26) 

where H(t) is the Hamiltonian (11) discussed in the pre- 
vious section. Fundamentally, the equation describes op- 
tical excitations acting as the external force driving the 
system out of equilibrium (the last three terms in (11)), 
and the dynamical couplings between QD and electrodes 
returning the two subsystems towards a state of mutual 
equilibrium (via both Hj< and V). Assuming that the 
electrode stays in equilibrium, we now reduce this equa- 
tion to the description of excitation, dephasing, and re- 
laxation of the QD alone. 

A. Reduced density matrix dynamics 

Recent experiments have pinpointed subtle FES ef- 
fects in optical spectra of quantum dots coupled to quan- 
tum wells 3 ' 41,42 . Thus it is crucial to include the ef- 
fects of FES in our theory. In order to achieve this, a 
special interaction picture must be constructed to fully 
account for the non-perturbative scattering of the elec- 
trode states in response to the electrostatic potential of 
the QD. Our approach is motivated by the one body 
formulation of the X-ray spectra by Mahan 28 , Nozieres 
and De Dominicis 43 , and the analysis of orthogonal- 
ity catastrophe by Anderson 31 . In the original FES 
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papers, the authors captured how the sudden shift in 
the Hamiltonian of electrons forming a Fermi sea af- 
fects the photoemission and absorption spectra in metal- 
lic systems 27 . Following this work, investigations into 
the FES in the photo-excitation of doped semiconductors 
were also carried out 24 . We also note the studies of the 
FES within the context of pump-probe experiments sen- 
sitive to the coherent non- linear optical response of doped 
semiconductors 44,45 . All these works focus on changes in 
the lineshape of the optical spectrum of a Fermi sea. 

Here we explore the consequences of these fundamental 
effects in the exciton dissociation across a QD-electrode 
interface. We formulate a master equation for an elec- 
trode coupled QD, and show how the effects of FES 
can be introduced in this theory, and how they can be 
re-captured naturally in the time-dependent couplings 
defining the resulting equation. Thus the FES becomes 
an integral part of the dynamical map that propagates 
the state of the QD towards equilibrium. The cou- 
plings in which the FES appears take the form of corre- 
lation functions bearing many similarities to the results 
of Nozieres and De Dominicis 43 . However, our equations 
address a very different physical scenario, and they are 
applied without making any simplifying assumptions on 
the spatial profiles of the different electrostatic poten- 
tials created by QD states. We also develop one-body 
formulas for calculating and interpreting these correla- 
tion functions. We will first discuss the derivation of the 
equation of motion for a general electrode, and then spe- 
cialize to the case of a Fcrmionic reservoir. 



In the Schrodinger picture, the reduced density matrix 
defined only over the QD states can be obtained by trac- 
ing over the electrode degrees of freedom. However, 
we begin with a slightly different version of this procedure 
and relate the reduced density matrix in the Schrodinger 
picture to the via the equation, 



P(t) = Tr L 



o -iH t/h 



me 



iHot/h 



(29) 



The FES arises when many body states are subjected 
to unitary rotation by two (or more) different Hamilto- 
nians. The exponentials in the above formula accom- 
plish this exactly, and allow us to expand the remain- 
ing interactions between QD and the electrode pertur- 
batively. From the chemical physics perspective, the ex- 
ponentials take into account the bath-induced random 
fluctuations of the QD energy levels, which cause the phe- 
nomenon of pure dephasing (decoherence without popu- 
lation relaxation) 48,49 . 

Returning to the general derivation, we write the 
Schrodinger equation within the interaction picture as 



#r (*),£(*) 
4 r) (*).£(*)' 



V c (t),i(t) (30) 



4 P) (*).£(*) 



The formal solution of this equation in the absence of an 
external electromagnetic field may be written as, 



1. Derivation of the general form 

From the full Hamiltonian defined in (11) in which the 
Coulomb interaction is given by (19), we take three con- 
tributions to construct the interaction picture by defining 
a Hamiltonian. 



Hn — Hj 



Hr 



(27) 



where n is a label for the QD states. The last term 
in the expression above scatters electrode states via a 
potential that is conditioned upon the QD state. This 
term is also the key to capturing the FES effects, but 
it complicates the interaction picture by yielding a non- 
perturbative coupling between the system and the bath. 

An analogue of this way of partitioning the Hamilto- 
nian has been used in the chemical physics literature in 
the past 46-48 . Such a reference system leads to the so 
called "modified Redfield" approaches 46,47 . The crucial 
difference here is that our bath is composed of fermionic 
excitations, which alter both the physics and the for- 
malism compared to the bosonic vibrational degrees of 
freedom at play in the previous work. 

To proceed with our analysis, we define the full density 
matrix within our interaction picture as, 



&t) = e iBot / h £(t)e- iHot / h . 



(28) 



i(t)=T+exp (-< / dt'Mt') U(0) 



(31) 



The symbol T + enforces time ordering such that 



T_|_ exp 



f dt'Ji(t') \ =1 -if dt'Ji{t') 



dt' 



dt' Jt^Mt") 



The superoperator ^Ji{t) in (31) acts on an arbitrary op- 
erator O as, 



Ji(t)0 



1 r 



H(t)-H ,O 



For our discussion below we also define superoperators 
for the radiative and Coulomb perturbations, 



Mt)o 



H T {i),0 



Jc(t)0 = - [V c (t),0 

To proceed, we assume that £(0) = p(0) <g> return to 
the Schrodinger picture and take the trace of (31) over 
the electrode states, 



p(t) = G(t,0)p(0), 



(32) 
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where G(t, 0) is a propagator for the reduced density ma- 
trix, 



G{ I A) ) = { e' lCot T+ exp (-i j dt'Ji(t'fj ^ . (33) 



We have introduced (-) L as the trace over electrodes, in- 
cluding £H. The superoperator e~ zCot is defined by its 
action on O as, 

g — iCatQ _ e — iH tQ e iH t 

To manage the subtleties of present choice of interac- 
tion picture, we explicitly find its matrix representation 
in the vector space of pairs of QD states by arranging the 
density matrix elements in a column vector with some ar- 
bitrary but fixed order. We define G as a matrix acting 
over vectors in this space with the matrix elements, 



(|o> <c|) 



t e -iC t 



T+ cxp ( -i J dt'Ji(t')\ 5R® |o'> (c'| 

We obtain the matrix elements G ac - a 'c' (t) by expanding 
the evolution operator on the right hand side to second 
order. Then, denoting its matrix form as G(t), we obtain 



G(t) 



= D(t) 



l- l D- 1 (t) f dt'(e- iCot J I (t')) 1 
Jo 



(34) 



— iCot 



dt' / dt" "JifflJiif) 



Here D(t) is defined to be a diagonal matrix over the 
same space as that of G, and its elements are given by, 



£W'&'(*) = Rab{t)e- iUabt 5 aa ,5 l 



w 



Rab(t) = Tr L [■ 



,iK b t-iK a t 



(35) 
(36) 



Here we have defined operators K a that act only on the 
electrode degrees of freedom, but are conditioned on the 
QD state, 



K„ 



1 



(H L + V a , 



(37) 



The general expression (34) for the propagator has also 
been derived earlier by Golosov et. al. is for a two- 
state system of electronic degrees of freedom coupled to 
bosonic nuclear motion. Below, we will specialize this 
propagator to a bath of fermions, and provide mathemat- 
ical details for the important distinctions with respect to 
a bosonic reservoir. 

To arrive at the equation of motion, we differentiate 
(32) with respect to t, set p(0) = G _1 (i, 0)p(t), and use 
(34) expanded to the order shown there. We thus obtain 



an initial value problem with the dynamics governed by 
a time-convolutionless master equation, which we write 
in the following form, 



d 

dt 



p(t) = M(t)p(t), 
P(0) = Po, 



(38) 



where p is now a vector obtained by re-arranging the 
matrix elements of p(t) as mentioned above. We express 
the time-dependent mapping in (38) as, 



M(t) = D^D- 1 ^) 



(39) 



+D(t)| {P(i) + B(t) + C(t) - V(<) j D- 1 ^). 

where P = Pc + Pt> and X signifies the derivative with 
respect to time. The form of this expression is motivated 
by the separate physical processes that are described by 
each of its terms, and the general expressions for these 
terms in the superoperator form are as follows (see (35) 
above for definition of D) 

V j {t) = -iD-\t)J t ti>(e^*J j (t , )) L , j = C,T (40) 

B(t) = -D- 1 (t) j dt' j dt" (e- lCat J T (t')J T {t")){Al) 

C(t) = -T>- 1 (t) j dt' j dt"(e- lCat J c {t')Jc{t")) L 

(42) 

Returning to (39), the first term of M(t) describes the 
decoherence caused by sudden switching of the QD po- 
tential, the general trends of which will be discussed in 
the following paragraphs. Note that T)(t) is a diagonal 
matrix by definition in (35). The second term in (39), 
arises from the first and second order contributions to 
the propagator, systematically included in (38) to sec- 
ond order. The first order contribution P(t) includes 
both hybridization and Coulomb interaction terms, al- 
though for the case of a Fermionic reservoir, only the 
latter will be non-zero. The second order contributions 
B(f) and C(t) respectively capture the hybridization and 
Coulomb interaction terms. 

2. Discussion and specialization to a Fermionic reservoir 



Let us now turn to the special interaction picture trans- 
formation defined by (28) and discuss its subtleties. We 
begin by taking the matrix elements of this equation be- 
tween two QD states, \a) and |6). Due to the fact that 



o iH t/h 



a iu a t+iK a t 



a), and similarly for \b), we find 



Pab(t) = Tr; 



I i(t)\b) 



iK,,l 



Therefore, unless is in the form of a product of the 
density matrices of the two sub-systems, a simple relation 
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does not exist between the reduced density matrices of 
the QD in the two pictures. In the derivation above, 
we have developed the dynamical equations under the 
assumption of a direct product form at t = 0. Therefore, 
it is instructive to analyze the consequences of Eq. (29) 
with the product form, i.e., 

f(t) = p(t) ® m, (43) 

where *H is any admissible density operator within the 
Hilbert space of the electrode states. The relationship 
between the two pictures is more complex than it is con- 
ventionally, 

Pab(t) = RabWPabWe-^-""*. 

Here, by p we mean the QD density matrix in the inter- 
action picture. Thus, in addition to the coherent oscil- 
latory factors arising from the QD states alone, there 
is an additional complex-valued multiplicative factor, 
Rab(t) ( see (36)), in the transformation from interac- 
tion to Schrodinger picture. The function R a b(t) is a 
manifestation of the Anderson orthogonality catastro- 
phe (AOC), which together with the Mahan exciton de- 
scribes the Fermi edge singularity effects 30 . To make 
the link with AOC more explicit, we pick a basis set 
{IV'a)} f° r the electrode states in which the operator 9\ 
is represented by a diagonal matrix with elements U\ a . 
Let | ipa(t)) = e lKat \ip a ), and write the function R a b(t) as 

iU(«)=5>«(^(*)|C(*)>- (44) 

a 

As pointed out by Anderson 31 , owing to the macroscopic 
size of the electrode, the overlap of the two states rapidly 
decays for a finite scattering of single particle states. We 
may also view the function R a b(t) as the average deco- 
herence caused by the electrode states, where the latter 
act effectively as a measurement distinguishing between 
the coherently superimposed states \a) and \b) of the QD. 

The solution to the equations of motion beyond the 
point at which R a b(t) vanishes, or crosses zero, can be an 
extremely poor approximation to the correct solution. To 
understand whether this presents a difficulty in our the- 
ory, we note that the leading contribution to AOC arises 
from differences in monopole moments of the initial and 
final potentials. Thus we expect the decoherence due 
to this mechanism to be very weak between QD states 
of the same net charge, so that R a b{t) ~ 1 within each 
class of states introduced above. When there is a charge 
transition, the AOC function decays as a power law at 
low temperature, or as an exponential at high tempera- 
tures. In cither case, the function does not vanish exactly 
within finite time except in the case of strong coupling 50 . 
To exclude strong coupling from the present scenario, we 
note that the coupling in our Hamiltonian is the Coulomb 
interaction, which is always screened by the electrode. 
For situations like those considered here the magnitude 
of the coupling is far from that of the strong coupling 
regime. In fact, we have verified this by explicit calcu- 
lations of R a b(t) for typical values of matrix elements in 



our Hamiltonian. Thus we proceed assuming that the 
functions R a b{t) decay but do not vanish exactly within 
the relevant temporal window of the dynamics. 

We now consider the consequences of specializing the 
above general formulation to the case of a Fermionic 
reservoir, which is the focus of the present work. We 
assume that the state of the fermionic bath representing 
the electrode is that of a normal metal, and described 
by a mixture of states $ a (N) where a is a state index 
for TV-particle many body states. Under this assumption, 
the annihilation operators entering H in (27) possess the 
property, 

c vh \* a (N)) = \&(N-l)), (45) 

for the single particle state \vk) having a finite occupa- 
tion in the many body state |$ a (iV)), and \&'(N - 1)) 
being (an un-normalized) N — 1 particle state. For such 
a state, the first order term in (34) vanishes whenever it 
corresponds to hybridization coupling. To see this, write 
the trace as a sum over the many-body bases \Q a {N)) 
and consider the matrix element between QD states \a) 
and | b), such that |6) has one extra electron relative to 
\a). Then the expectation value, (e~ lCot J'i(t')), in the 
first order term consists of two terms. One of these is 
proportional to the following sum, 

^SH Q (7V) ($ Q (N)\ e -« f »(*-* , ) c t fcC **»C*-0 |* (JV)> , 

and the other is obtained by c^, fc — > c„fc. In the equa- 
tion above, 9\ a (N) is the probability of state & a (N) in a 
grand canonical ensemble. Since, the operators K a and 
Kb do not change the total number of particles, the result 
of the above formula is an overlap between A^-particle and 
N + 1-particle Fock states. Therefore it vanishes, and so 
do all odd order terms in (34). 

In the language of our formalism, this implies that 
Pt = 0. We remark that this is a consequence of (45), 
and the fact that hybridization involves an odd number of 
creation annihilation operators for electrons. If the elec- 
trode were, for example, in a BCS state then the previous 
sum over states would also involve states in coherent su- 
perposition of different N, resulting in a non-vanishing 
expectation value of c* uk in general. 

Furthermore, when considering the contribution of 
Coulomb coupling via J^Ci the first order term, Pc in 
(34) does not vanish. However, since Coulomb coupling 
does not change the charge state of the QD, this term 
can be understood as a Hartree energy correcting for the 
fact that electrode states are defined in the presence of 
a neutral QD. Thus for a Fermionic reservoir, Pc(t) in 
(39) is the off-diagonal Coulomb potential matrix. It is 
discussed below in Sec. Ill C with its definition given by 
(59). We also note that this term does not affect popu- 
lations, but describes only the dynamical re-organization 
energy within the electrode during coherent oscillations 
between different charge states of the QD. 

For an electrode at equilibrium, the charge and energy 
transfer processes may couple at third order. This cou- 
pling would affect the rates of charge transfer such that 
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an excited electron or a hole in a state with low escape 
rate may exchange energy with the electrode and jump 
to a state with larger escape rate. This modification of 
the charge lifetime of the QD may be expected for closely 
lying hole levels. However, since the Coulomb interaction 
is small due to screening, and tunneling is exponentially 
suppressed by increase in junction width, this regime may 
be an exception rather than a rule. Thus we do not pur- 
sue it in the analysis below. 

We now stay within the confines of a Fermionic reser- 
voir, construct the expressions for B, Pj, and C matri- 
ces, and discuss their physics. From the physical insight 
gained into these matrices, we are also able to obtain 
useful approximations that simplify their numerical eval- 
uation. 



Figure 2: Double sided Feynman diagrams showing the 
effect of the three terms in the dynamical hybridization 
couplings. The red arrows represent the addition and 
removal of electrons from the QD respectively, and the 
two vertical lines represent the "ket" on the left and 
"bra"on the right. 



B. Charge transfer rates 

We first write the matrix elements of B(t) in terms 
of the hybridization operators to facilitate the connec- 
tion with electrode correlations, and then evaluate the 
matrix elements as shown below. The matrix elements 
B aC ;a'c' (t) follow from the general expression (41) special- 
ized to a Fermionic reservoir, in which the term = 
as discussed above. Thus we obtain, 



B ac .y c ,{t) (46) 

^2{tcc> (T ca (t) T ah {tx)T* a ' '(ta)3t\ e l{u - btl+UJ ^' t2) 

b 

+Saa> (r c ' fc (t 2 )r bc (ti) r ca {t)mj e ^'^+^ ti ^ 



h 2 



dti / dt2 



T c ' c (t 2 ) T ca (t) T aa '(h)m) e < K<»'*i+"c'ct2) 



where uj ac = u a — uj c . The matrix element B ac . a > c >(t) 
describes scattering from the "initial state" p a i c i to the 
"final state" p ac . Here the phase factors due to the QD 
energy levels are shown explicitly and the time dependent 
hybridization operators reflect the action of K a in (37). 
In addition, we have introduced what we call the pure 
dephasing operator, 



e iK c t e -iK a t 



(47) 



which does not generate any charge transfer, and only 
affects coherences by accounting for the FES effects for 
oscillations between states of the QD with different effec- 
tive Coulomb potentials. We remark that it is straight- 
forward to verify that the matrix elements, B ac;a / c /(t) 
obey the sum rule, 



B,, 



t (i) = - B n 



(48) 



which is quite general, and in turn ensures that the sum 
of all the rates for the population to relax from a state 
\m) equals the total decay rate of this state. 

Returning to (46) , we note that since there are no hy- 
bridization operators of the form T cc (t) in the entire 
Hamiltonian, a given density matrix element is acted 
upon by either the first two terms of that equation or 
the third but not both. However, the three terms are not 
independent, and satisfy sum rules due to the conserva- 
tion of total particle number by the underlying Hamil- 
tonian. These terms may also be interpreted as general- 
ized scattering "in" and "out" rates for populations and 
coherences. We depict the effect of these rates on the 
density matrix in Fig. 2 using double sided Feynman di- 
agrams. As shown there graphically, the first two terms 
couple only coherences to populations, while the third 
term provides an additional pathway for coupling popu- 
lations that differ by one electron. The same effect would 
occur at a higher order in the form of the third diagram. 
Note that 7o(t) does not have a representation in terms 
of these double sided Feynman diagrams because it does 
not change the state of the QD. 

From the definition (15)-( 18) of the hybridization op- 
erators, and for an electrode consisting of an electron 
reservoir in equilibrium, we find 



Bac;a'c' (^) 



(49) 



h 2 



dt' 



dco 
2tt 



7c'c;aa'(w, t' ] t)S{uj - UJ a , a ,t') 
+7a>a;ce> (w, t'; t)S(u> ~ Ujcrf) 

Sa'a X cb;bc> ( W > t) S ^ ~ ^cb, t') 
b 

b 

The derivation of this formula is provided in Appendix B, 
and we have defined here a function. 



sin(wt) 2zsin [wj] 



+ 
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which accounts for the initial condition defined at a fi- 
nite time, and allows us to work with Fourier transforms 
with respect to the initial time at — oo (see Appendix 
(B)). The functions 7 c < c;oa /(u; ,t';t) and X c a a b . ba , (ui, t'; t) 
are a generalization of the spectral functions and defined 
as Fourier transforms of causal response functions, 



+00 



lc'c-aa' (w, t';t) = -21m / dt"T dc . aa> (t" , t'; t)e iut " , 



-\-oo 



\ c a U a ,(uJ,t';t) = -21m / dt"AZ :ba ,(t",t';t)e^' 



so that the Fourier transform integrals can extend over 
the entire real axis. The integrands in the previous ex- 
pressions are the correlation functions appearing in (49) ; 
the superscript in A identifies the indices on the pure de- 
phasing operator (47), while these indices are implied by 
subscripts in the function T. The function A captures the 
first two correlation functions in (46), while T captures 



the last one, and they are defined as, 

Kl% ba ,(t",t';t) = -i@{t' - 1") (r Q ca (t)T ab (t')T ba '(t")m) , 
r cW (t",*';i) = -i®(t' - 1") (r c ' c (t')T ca (t)T aa '(t")m) . 

(50) 

In contrast to conventional correlation functions for a 
single particle propagation, these two correlation func- 
tions contain three time arguments. Their dependence 
on the third argument arises from the operator 7~Q a (t), 
which differs from unity only when a =/= c in general. 
Furthermore, physically, it is only significantly different 
from unity when the difference between the potentials of 
the two states is large enough to cause significant AOC. 
Thus the third argument describes the shakeup of the fi- 
nal electrode states when the QD oscillates between two 
different charge states. This is fundamentally different 
than the processes described by the first two arguments. 
Specifically, the time difference t" — t' relates to the parti- 
cle absorption/emission spectrum of the electrode in the 
presence of a QD. The average time (t" + t')/2 relates to 
the memory of the initial state potential of the QD. In 
our notation, we use a semi-colon to set apart the two 
different kinds of time arguments. 

The full mathematical analysis of correlation functions 
Eq. (50) is outlined in Appendix B. To gain physical in- 
sight into the results therein, we focus on the process of 
electron transfer from the QD to the electrodes. Then 
the pertinent correlation function is denoted by the su- 
perscript ">" and defined as, 



r c % aa< (iV;*) = ^e(*'-i") £ T%T™ k , ( e ^*' c ^ e -^(t'-t) e ^( t »- t ) c t /fe/e -^ a , t » e -^^ ; (51) 

vv'kk 1 



where we have made the physical assumption that *H = 
exp(— (3Hi,)/Z with the usual partition function Z = 
(exp(-0H L )). 

The correlation function under the sum on the right 
hand side of (51) may be interpreted as a thermal av- 
erage of overlaps between two electrode states evolving 
under different potentials. To see this, take a many body 
electrode state with N electrons, |<fr(JV)), which would 
appear as a ket when computing the expectation val- 
ues in (51). First consider the case where 7^, ca (i) ~ 1, 
which occurs when the states |a) and |c) have the same 
charge. Then all the correlation functions above reduce 
to functions with two essential time arguments. The state 
|$(A)) initially evolves under the potential V a ' a > until 
time t" , and at this time an electron is injected from the 
QD into a single particle state of the electrode with in- 
plane momentum k' . Within this scenario, the QD can 
either be in a negatively charged state or in an exciton 
state immediately before t": it cannot be in the ground 



I 

state. After tunneling, the QD potential switches to V aa 
and the state is evolved back to the initial time (t — 0). 
Let us call this state |<I>"(A + 1)) to denote the fact that 
the state corresponds to time t". Similarly, the bra form 
of |$(A)) evolves under V c i c i until time t', when an elec- 
tron is added to it in the single particle level k, and the 
resulting N + 1 electron state is evolved back to the ini- 
tial time under the influence of the QD potential V cc . We 
label the resulting state \&(N + 1)). 

With these definitions, we see that the corre- 
lation function in (52) is equal to the overlap 
($'(A + 1) I $"(A + 1)), the trace in that expression be- 
ing a thermal average over all the many body states, and 
the summation over k and k' is a summation over all 
possible in-plane momenta into which the electron can 
be injected. The weights for the latter summation are 
the tunneling amplitudes, which give the probability of 
such an injection to occur. We may interpret in the 
same way the correlation functions in which c^, is to 
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the left of Cfc- As they describe the evolution with an 
electron removed from the initial state, the overlap is 
(&(N- 1) | $"(JV-1)). 

It is now clear that if V c > c ' an d V a ' a ' are weak (for ex- 
ample, if their source has a vanishing charge), the prop- 
agation of |$(AT)) until the time of adding the electron 
is identical in both |$'(JV + 1)) and |$"(JV + 1)). There- 
fore, the origin of time loses significance and the resulting 
correlation becomes a function only of the time difference 
t' — t" . It thus reduces to a correlation function for a sys- 
tem in equilibrium or in a steady state, where the effects 
of initial conditions have decayed. Therefore, the depen- 
dence of the generalized spectral functions on t' in (49) is 
in proportion to the memory of the QD potential in the 
initial state. Similarly, it follows from (50) that the de- 
pendence on t is in proportion to the difference between 
potentials of the superimposed final states at time t. 

It follows from above that if only the monopole contri- 
bution is significant, the correlation function would de- 
pend only on the time difference, t' — t" so long as the 
initial and final state of the QD have the same charge. 
The Fourier transform of the correlation function with re- 
spect to this time difference then defines a conventional 
spectral function describing single particle absorption or 
emission as a function of energy. When coupling to the 
dipole and higher order potentials is strong enough, co- 
herent oscillation between the states of different multi- 
poles would, in general, lead to this spectral function 
evolving as a function of the average time (t! + t")/2, 
thus exhibiting non-Markovian behavior. 

We now summarize the mathematical expressions that 
can be evaluated with our physical model. In particular, 
we specialize to the case of quadratic Hamiltonians, and 
follow Blankenbecler et al 51 and Hirsch 52 to develop ex- 
plicit expressions for these correlation functions in terms 
of a matrix inverse and matrix multiplication in the single 
particle basis for the electrode states. Thus for the corre- 
lation function in which an electron tunnels from the QD 
to the electrode, the analysis presented in Appendix B 
yields, 

r$ c;aa , (t",t';t)= (52) 

-ie(t' - t")r c , c , aa , (*", a, t) T^Tti 

v.v'kk' 

(uk\ N c , c , (-*', 0) [i + Eg,,, (f, f"; t)N c , c , (-f , 0) 
v'k>), 



e -iK a ,t' e iK a ,t 



and 

>ca ( ±11 if 



A^(f'^)= ( 53 ) 
-iQ{t' -t")r ca ., w {t,t\i!) £ Tf k T% 



'kk' 



uk\ N ac (t - t', -t) [i + Eg, (f, t"; t)N ac (t - t', -t) 
u'k') . 



The superscript ">" in these definitions indicates that 
the correlation functions correspond to the propagation 
of a state in which an electron is added to the electrode 
(see Appendix B). The pairs of indices on the correla- 
tion function are constrained by the type of correlation 
function, "lesser" and "greater", in accordance with (B5). 

In these expressions, the evaluation of the correlation 
functions has been divided into two parts. First, we have 
defined a more general function to capture the AOC ef- 
fects for four potentials, 



^c' c;aa' (t 3 ^ ) ^) 



(54) 



Z 



e iK c ,t' e -iK e (t'-t) e iK a (t"-t) e -iK a ,t" e ~0Hi 



Second, the part that explicitly depends on the hybridiza- 
tion functions is now reduced to the proper combination 
of matrices that represent each indicated operator in the 
natural single particle basis for the electrodes, \vk) . The 
matrix N QC (i,r) is defined to be, 



-iK a t -iK c r n j_ „-mLY- i - e iK c r e iK a t 



N ac {t,T) = e- lK ^e- lK ^ 1 + e" 



. (55) 

In the \ vk) basis, the central term is simply the statistical 
weight for empty states 



^vk,i>'k' 



l + exp{-(E uk -(x)/k B Ty 



(56) 



e -iK e ,t' e iK a ,t' 



Thus, the process for electron transfer into the elec- 
trode is proportional to the empty states available. How- 
ever, the N QC (i, t) includes the time-dependent shake- 
up processes induced by the potential of the QD states 
I a) and |c). In turn, this influences the distribution of 
states available to receive the electron, opening addi- 
tional states below or blocking states above the Fermi 
level. 

The matrices S^, j = 1,2, in (52) and (53) are com- 
mon to both emission and absorption, and are defined 

as, 

S^ ; l,(i',i";*) = e^ iKc ' t ' e lKa,t " e~ iKat " To b ' (t)e iKct ' — I, 

(57) 

4 2 l(t' 1 t";t) = e- lK ^'% bb ' (t)e lK ^' t " e-' lK ^" e lK ^' - I. 

(58) 

For the correlation corresponding to an electron tunnel- 
ing from the electrode to the QD, expressions similar to 
(52) and (53) exist (se_e Eqs. BIO and B12). In those 
correlation functions, N is replaced by the compliment, 
N = I - N. 



C. Energy transfer rates 

We now turn to the expression for the matrices P(t) 
and C(<), which are the first and the second order terms, 
respectively, that enter in the cumulant expansion of the 
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propagator (34). The first order term deriving from the 
Coulomb interaction, Eq. (40), and specialized to the 
Fermionic reservoir, reads 



Pac;a>c>{t) = -iD-}{t)- 



-iuj ac t 



h jo 

iu) nn ft' 



Sec (T ca (t)V aa ,(t')K)e 
5aa-(Vcc{t')% ca {t)V\) 



(59) 



Recall that the diagonal terms (e.g. V aa ) do not ap- 
pear in the Vc operator, having been put in Hq. As 
shown in (B13) the above expression for P ac -a'c' may be 
expressed in terms of a Hartree energy matrix. 



A aa ,(r) = (l-S aa ,)(V aa ^ aa '\T)) 



(60) 



the time-dependence of which arises from the evolution 
of the electrode state under the influence of the average 
of the potentials of the QD states \a) and \a'), which we 
write as 



£K (a<,,) (t) 



me 



By writing the time-dependence in this manner, we fully 
capture the monopole contributions to pure dephasing 
(see Appendix B 2) . The effects of first and higher mo- 
ments of this density can be ignored due to the Friedel 
sum rule 29 , and have been verified by our calculations to 
be small. We remark that this approximation is a direct 
consequence of the charge conservation implied by the 
off-diagonal components of the Coulomb matrix in (19). 
Furthermore, the Hartree correction arises only in the 
coupling of coherences to each other and to populations 
due to the fact that it is of first order and that the diago- 
nal elements of its underlying Coulomb potential matrix 
have been removed. 

We now turn to the second order term, denoted by the 
matrix C(i), which generates energy transfer between the 
electrode and the QD. The full expression for C(t) is iden- 
tical in form to that for B, except for the appearance of 
the Coulomb interaction operators V ac (t) instead of the 
hybridization operators Tacit)- However, instead of start- 
ing with the full form, we derive the expression for C by 
neglecting the effects of AOC altogether 53 . These effects 
can occur only at third order in the Coulomb interac- 
tion between the QD and the electrode. Furthermore, 
retaining these effects amounts to calculating the dielec- 
tric function of an electrode driven out of equilibrium by 
the fluctuations in the QD potential. We expect this to 
be a very weak effect due to the screening of this poten- 
tial, and the macroscopic size of the electrode. Further- 
more, since the Coulomb interaction does not change the 
charge of the QD states, changes in the potential induced 
by Coulomb processes are much weaker than hybridiza- 
tion. The neglect of AOC in energy transfer simplifies the 
expression for Pc(t) by reducing it to Pc?(t) — (Jcif)) l- 
This also allows us to write the combination C — P 2 /2 



in (39) in terms of the fluctuation of potentials around 
their average values, with (42) giving the general expres- 
sion for C(t) 

Thus for states \a) and \a') carrying an equal charge, 
and \d) and |c) carrying an equal but not necessarily the 
same charge as \a), we obtain, 



a 



ac;a c 

i 



it) 

■iu) ac t 



(61) 



t rti 

dti / dt2 
Jo 



E 



5 dc (AK^OAVWfe)^) e^^+^'t*) 
' AV c , b (t 2 )AV bc it 1 )m) e^eti+wc'^)' 



1 



t rt 

dti I dti 
li Jo 



A^ c (t 2 )AKa'(ti)5H)e 



The three terms in this expression act on the density 
matrix in a way similar to the corresponding three terms 
of the tunneling process. Here we have defined the op- 
erator for the fluctuation of V a b around its mean value 

as, 

AV ab (t) = e iH L t/h^ abe -iH L t/H _ (y^ 

Note that we have neglected the effects of AOC due to the 
QD potential already in this definition; the subscripts on 
these operators only identify the matrix elements when 
(20) is substituted for performing calculations. We re- 
mark that the matrix elements C ac;a > c >(t) also obey the 
sum rule analogous to (48) 



mm;mm 



it) 



Following the derivation of 7, 
now define the functions 



Xc 



(*)• (62) 
(uj,t';t) above, we 

2Im(l-^ c ) (!-*„/„) J] (63) 



V c :, v - c y{-q)X, 



iq,^)V a ,^a'^'iq), 



where X vv i- tlJll j,i{q,U3) is the renormalized density-density 
correlation function, which is generalized to include inter- 
subband transitions. By renormalization, we mean that 
a summation over an infinite series of particle-hole exci- 
tations is performed between the times of the two interac- 
tions. Thus if X® v ,. , (g, lu) is the response function of a 
non-interacting gas, and corresponds to a single particle- 
hole excitation, then 



>(q,w) = X°„, ,(q,u) 



X 



, iq, uj)V aa/ -pp> (q)Xpp>. aa > (q, u), 

(64) 

where V aa '-BB'iq) is the effective Coulomb interaction 
describing the momentum exchange q between particles 
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scattering from bands a', j3' and into the bands a, (3. This 
matrix is an effective interaction because it accounts for 
the static screening by the bulk substrate beneath the 
surface. When this substrate is a metal, the large plas- 
mon frequency ensures that the bulk response may be 
considered instantaneous, which creates a static screen- 
ing of the Coulomb potential of a surface electron, and 
the result defines the two-particle interaction for the sur- 
face modes. Similarly, the effective bulk dielectric func- 
tion is also static for insulators or wide-gap semiconduc- 
tors. The dynamical screening by the bulk may be im- 
portant only for narrow-gap materials. 

Following the mathematical steps outlined in Ap- 
pendix B we obtain, 

:C ;rt '(t) (65) 



dt 
1 



e ~iui aa t e i(uj a 



0/ 



2h? 

Xc'c;aa' 
~\~\a' a\cc' 



d/jj 



Xc'b;bc{u)S(uJ - OJbc',t) 



b 



Xa'b;ba(u)S(u - U>ba',t). 



This expression forms the basis of our calculations of 
energy transfer driven by Coulomb interaction between 
the QD and the electrode. We remark that the elec- 
tron susceptibility in this expression may be calculated to 
any level of sophistication within the computational con- 
straints. Note that we did not include electron-electron 
interaction within the electrode explicitly in the Hamil- 
tonian (11). These interactions enter our theory via the 
correlation functions, which depend on elementary ex- 
citations that transfer energy between the electrode and 
the QD. Thus they can instead be taken into account fully 
by constructing appropriate dynamical dielectric func- 
tions and single particle Green functions. 

The analysis presented so far applies only to the intrin- 
sic couplings in the system, which lead the sub-systems 
to mutual equilibrium. The radiative interactions that 
drive these sub-systems out of equilibrium can be con- 
sidered with the additional light-matter coupling shown 
in (23)-(25). The full development of optical response of 
these systems can be developed based on the above for- 
malism. By computing the linear and non-linear optical 
response of surface coupled quantum dots in this way, 
one can develop ways in which advanced spectroscopic 
techniques can yield experimental measurements of the 
sub-system couplings. However, this is beyond the scope 
of the present paper. 



IV. APPLICATION 

We now apply our model to the system shown schemat- 
ically in Fig. 1 in the beginning of the paper. The goal 



of this section is to illustrate the calculations of the main 
quantitites in the formalism in the order it has been pre- 
sented above. We first calculate the microscopic Hamil- 
tonian from a model of InAs QD, and an ab-initio model 
of a Au[lll] electrode. We use the electronic states com- 
puted from these models to obtain couplings between the 
subsystems. We then illustrate the use of the results of 
these calculations as input to our dynamical theory. In 
particular, we calculate the charge and energy transfer 
rates, and the non-radiative exciton recombination rates. 
In these calculations, we show the impact of FES and im- 
age effects on the dynamical charge and energy transfer 
couplings. 



A. Quantum Dot states 

The QD is modeled as a dielectric sphere of InAs with 
radius, a = 2.0 nm, and dielectric constant €qo = 15.0. 
We set the distance between the QD center and the im- 
age plane of the electrode as h/2 — 2.625 nm. Following 
Chulkovef al. li we set the location of the image plane 
plane to be 2 A above the top atomic layer of the elec- 
trode. To model the QD states, we use the effective mass 
approximation, in which we include only the lowest con- 
duction and heavy-hole bands of bulk InAs. Thus we 
use (7) in the form 



2mo m j( r ) 



V + V{r) + S(r) - e t , 



Vn{r) = 0, 



_ (66) 

where rrij is piecewise continuous with rrij = 1 in the 
space between the QD and the electrode. Inside the ge- 
ometric boundary of the QD, it is equal to the effec- 
tive mass of the conduction band or a heavy-hole band 
for electron (j = e), or hole (j — h) states respec- 
tively. Inside the QD, we set m e = 0.03, and account 
for the anisotropy of the InAs heavy hole bands by set- 
ting rrih = 0.52 parallel to the electrode surface and 
rn/i = 0.33 perpendicular to it. Note that this type of 
modeling implicitly sets the smallest spatial scale to be 
the lattice spacing of the QD such that its boundary is a 
shell of zero thickness. 

The potential V{r) is a square well potential of a spher- 
ical QD with V = outside the QD and V = -5.0 
eV inside. The value inside is typical of the workfunc- 
tions of semiconductor materials, and in our discussion 
of dynamics below, we explore how the charge transfer 
rates depend on the alignment of the square well po- 
tential with the Fermi level of the electrode. The self 
energy, S, describes the electrostatic reaction field due 
to the polarization of the electrode and the QD surfaces, 
thus taking into account the image effects of both sur- 
faces. An exact calculation of this electrostatic potential 
is used, the details of which can be found in our recent 
publication 34 . Since our smallest spatial scale is larger 
than the lattice constant, we interpolate the electrostatic 
image potential of the QD across this width. Another 
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essential assumption of this electrostatics calculation is 
the uniform dielectric constant inside the QD, which is 
well-justified by several ab-initio calculations published 
in recent years 54 ' 55 . 

Using cylindrical symmetry, we reduce the problem 
to two spatial dimensions, normal and parallel to the 
electrode surface, and solve (66) numerically over a two- 
dimensional grid of (r, z) points using the method of finite 
differences. We use the the so-called ghost fluid method 56 
to subsume the mass discontinuity 34 . 

Our calculated total potential for the electron in the 
geometry described above is shown in Fig. 3. We see 
from this figure that the image attraction by both the 
QD and the electrode lowers the potential significantly 
in the narrow tunnel function around the surface normal 
passing through the QD center. This plays an important 
role in increasing the rate of electron tunneling. On the 
other hand tunneling is affected little for the states whose 
symmetry places a node in the wavefunction within this 
junction. 

We calculate the exciton state by employing our exact 
solution for the electrostatic polarization to determine 
the matrix elements of the effective two-particle interac- 
tion potential, V(r e ; r^), in which a full account is taken 
of the interaction of a hole with the volume density of 
the electron, and the surface charge it induces 34 . We 
then expand V(r e ;rh) in the basis of product states ly- 
ing below an energy cutoff, and increase this cutoff till 
the binding energy of the lowest exciton state converges 
to within 1 meV. The effect of electron-hole correlation 
on the resulting state is subtle but important. While the 
lowest exciton wavefunction is dominated by the product 
of the lowest electron and hole states, the mixing of states 
introduces corrections that can have noticeable effects on 
decay rates, as discussed below. 

The right hand side panel in Fig. 3 shows the en- 
ergy levels obtained after the exact diagonalization just 
described. The exciton level indicated on the figure is 
shown relative to the hole level of the QD . One may view 
this as representing the electron level correlated to a hole 
in the top valence level. 



B. Surface states of the electrode 

To compute the electrode wavefunctions at the sur- 
face, we first perform an ab-initio calculation (using the 
SIESTA program 57 ) for bulk Au, followed by a supercell 
calculation with 24 mono-layers of Au oriented along the 
[111] direction separated by the equivalent of 36 layers of 
vacuum. The surface unit cell is lxl, with no reconstruc- 
tion. The TRANSIESTA 58 code is then used to generate 
the Green functions for the surface Au layers via recur- 
sive decoupling 59 ' 60 from the bulk layers. From these 
Green functions, we determine the projected density of 
states (PDOS) in which the Shockley band 61 is identified 
as shown in Fig. 4(a). 

From the localized basis orbitals generated by SIESTA, 




Figure 3: Left: Potential energy of an electron in the 
equatorial plane outside the electrode. Right: Energy 
levels of electron, hole and exciton states. The exciton 
levels is placed such that the hole of the exciton resides 
in the top hole level. The dashed dotted lines refer to 
conduction (top) and valence (bottom) band edges. The 
thick dashed line is the Fermi level of the electrode. 



and the eigenvectors at the poles of the surface Green 
function, we calculate the wavefunctions for the surface 
states along this band. The eigenvectors are taken along 
the peak of the PDOS of the surface band. Fig. 4(b) shows 
the state at the bottom of the band superimposed on the 
atomic layer positions of Au[lll]. Since the dispersion 
of this state is predominantly parabolic, its exponential 
decay remains essentially constant along the band. This 
is because the rise in total energy of a parabolic band 
is cancelled by the rise in its kinetic energy parallel to 
the surface, thereby leaving the wavenumber along the 
surface normal fixed to its value at the bottom of the 
band. 



C. Couplings 

From these wavefunctions, and those of the QD, we 
determine the hybridization matrix elements Ht from a 
direct application of (A6). To calculate the Coulomb in- 
teraction matrix elements, we employ the procedure de- 
scribed by Eqs.(All-A14). The main input to this pro- 
cedure is the effective potential in which corrections due 
to the surface dielectric function of the electrode are ap- 
plied. In order to do that, we follow Pitarke et al. 62 
and divide the surface into two linear response systems. 
The first is the semi-infinite bulk, and the second is the 
two-dimensional electron gas formed within the Shockley 
band. We let Xb an d Xs be the susceptibilities of the two 
systems. 

We define e^ 1 to be the dielectric function of the bulk, 
and follow Newns' work 63 in calculating it. Newns' cal- 
culation is based on the random phase approximation 
(RPA) within a jellium model of the electron gas, and 
the potential we calculate thus represents screening by 
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Figure 4: (a) Density of states of Au[lll] surface close 
to the r point (a), (b) The averaged potential along the 
surface normal. The superimposed red curve is the 
Shockley state. The atomic planes are indicated as 
dashed lines, while the image plane is shown as the solid 

line. 



given by 



W(q,z) 



V (q,z) - V (q,z) 



1+^(9) 



Q(z - z p ) 



+ V ( q ,z p ) 2£ ~y> Z Zp) 8(z p - z), (67) 

where e^ 1 (q, z — z p ) is obtained by multiplying the pre- 
factor in formula (62) of Newns' paper 63 by —q, and z p is 
the coordinate of the receded plane. When Vo(q, z) corre- 
sponds to the potential of the charged QD, W(q, z) enters 
the Hamiltonian matrix via (20) and accounts for the sub- 
strate response. The response of the surface gas appears 
explicitly in the dynamical rate expressions, and this 
asymmetry in our treatment arises because the dynamics 
depends only on frequencies that are much smaller than 
the bulk plasmon frequency. This allows us to use only 
the static response of the substrate, but since the surface 
acoustic plasmon branch extends to zero frequency, the 
full dynamical response of the 2D surface gas must be 
included. Furthermore, the coupling of the QD poten- 
tial is also much larger to the electron hole excitations 
in the surface gas than it is in the bulk, since the poten- 
tial is screened fully within a few atomic layers below the 
surface. The dominant coupling to surface states is also 
verified by the ab-initio calculation of electrode states as 
described above. 

Thus, having made the choice to let W(q, z) represent 
the entire bulk response, we now turn to the susceptibility 
Xs of the 2D surface electron gas. In the present imple- 
mentation approximating the general picture developed 
in Sec. IIIC, the energy transfer rates originate from the 
response of the 2D electron gas. We calculate the non- 
interacting x° using Newns' approach 63 , and use RPA to 
construct the interacting response Xs as in (64). In the 
RPA calculation, the electron-electron interaction within 
the 2D gas is screened by the substrate, but since the 
screening is only partial, collective excitations in the 2D 
electron gas can exist in the form of surface acoustic plas- 
mons 62 . We find the effective screened interaction using 
Newns' work, 



the semi-infinite bulk only. Due to the large plasmon fre- 
quency of the bulk, we replace £& by its static limit. We 
include the effects of the surface states based on the work 
by Pitarke et al e2 and Silkin et al 64 . In their simplified 
model, the surface states comprise a 2D electron gas lying 
in a plane a small distance \z p \ above the interface of the 
semi-infinite jellium representing the bulk substrate 62 . 
This separation ensures charge neutrality in the interior 
of the jellium medium, 63 and in the present case we place 
the 2D plane at z = 0, and set z p = —"in/Skp s» 2A 64 . 

Given Vo(q, z) as the planar Fourier transform of a bare 
potential, the potential screened by the jellium plane is 



v cS (q) = v(q) - v(q)e- 2 ^ 



1 



(68) 



Here the first term is the Fourier transform of the bare 
interaction within the 2D plane representing the surface 
electrons, and the second term is the interaction with the 
image charge in the substrate located in the plane a dis- 
tance 2\z p \ below the surface electron plane. Substituting 
v c g(q) into the RPA summation for the 2D response, we 
obtain the screened susceptibility, Xs(q,u), as 



X°(9,w) 



1 - v cS (q)x° s (q,uj)' 



In the left panel of Fig. 5, we have plotted a two- 
dimensional color map of the function Imx s (g, uS) in 
which the bulk screening is neglected. We have verified 
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that the acoustic plasmon branch in the a bsence of bulk 
substrate follows the expected trend, yj e 2 nq/ (2Meo), 
which holds at low energies for a 2D gas of quasi-electrons 
of mass M and density n. By including the screening due 
to the bulk substrate, we obtain the plot shown in the 
right panel. The plasmon branch in the result is signifi- 
cantly different as it has now acquired a linear profile that 
straddles along the incoherent pair continuum. Thus the 
bulk substrate shifts the spectral weight of the plasmon 
branch to a lower frequency and spreads it into the pair 
continuum. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 ! 
q [1/nm] q [1/nm] 



Figure 5: log 10 [— Im%] in the absence of the 
semi-infinite bulk (left), and in the presence of static 
screening by the bulk (right). Note the qualitative 
change in the dispersion of the acoustic plasmon 
branch, which in the presence of the bulk lies close to 

the incoherent electron-hole continuum, the upper 
boundary of which is marked by the thin dashed line. 
Thus the spectral weight of the plasmon contributes to 
the pair continuum. 



Let us briefly comment on when this shift can become 
important. Consider the processes in which a state n 
relaxes to a lower energy state m via the Coulomb in- 
teraction, in which the coupling potential is formed by 
the transition density <p* m (r)ip n (r). From (63) and en- 
ergy conservation, we expect that the region in the {q 1 E) 
plane that couples to this process must lie within the 
Fourier-Bessel transform of this potential along the q 
axis, and in the range of energy differences E n — E m 
along the E axis. Wavefunctions inside a QD of radius 2 
nm would generally yield the peak of the Fourier-Bessel 
transform to be around q ps 0.5nm _1 . The plot in Fig. 5 
then implies that E n — E m must be 0.4-0.5 eV apart for 
the states n and m to experience qualitative change in 
their coupled dynamics. As shown by our calculations 
below, this energy difference is much higher than the in- 
verse rates implied by the energy transfer matrix C(t) 
when the bulk substrate included. On the other hand, in 
the absence of the substrate, it is possible to approach 
the regime where these changes may significantly modify 
the time-dependence of energy transport between the QD 
and the electrode. We now turn to the effects of charge 
and energy transfer on the exciton and hole populations 
with the initial state in which the QD is prepared in the 
exciton state by photo-excitation. 



D. Exciton dissociation and Fermi edge singularity 

In this subsection, we discuss the charge transfer of 
the electron to the electrode, which leaves a positively 
charged QD containing one hole. Since the charge state 
of the QD changes from neutral to positive, we expect 
to see the effects of the FES in the tunneling rate. With 
reference to the full rate expression in (39), we focus on 
Mhh-xx so Dhh and D xx enter, both of which are simply 
unity, and the rate is simply, 



b x {t) 



h 



A brief description of the calculation strategy is as fol- 
lows. From (49), we obtain 



B 



hh:xx 



(*) = 



2 



duj 
2^ 



~fxh;hx{u)S(u) - Ul xh ,t). (69) 



The controlling energy scale is the excess energy of the 
tunneling electron relative to the electrode Fermi energy, 
expressed by u x h- We evaluate the above expression us- 
ing (52), (54)-(58). Since the dipole field of an exciton is 
weak and it does not affect the edge singularity by the 
Friedel sum rule 29 , we set [K x ,Hl] ss 0, which also sets 
N xx (t, 0) fa N. The definition of J x h;hx(w, f; t) and this 
approximation together imply that it is independent of 
the time arguments t' , t, which we reflect in the equation 
above by omitting the time arguments (see also the dis- 
cussion at the end of Sec. IIIB). Taking only the surface 
states as the final states in the charge transfer, we obtain 



l xh -M{u) = -iY. T $ T sV 



(70) 



fefe' 



+ 00 



dre lUT r xh , hx (T) x 



sfelN 



I + S^(r)N 



Sk') 



where s denotes the surface states, and 
1 

~Z 



^xh-hxV) 



r xh . M (r) = ^(e iKhT e- iK ^e-^), 



Note that (70) is different from the Fermi golden rule 
result in that it includes the effects of sudden switch- 
ing of potential to all orders in perturbation theory, and 
therefore includes the effects of FES exactly. This re- 
sults in dressing of all operators by the Coulomb poten- 
tial of the QD. The above result is a perturbation expan- 
sion with the dressed hybridization as the small parame- 
ter. Eliminating the FES effects is equivalent to setting 

r x h-hx( T ) = !) and ^xh°hx( T ) = °- Substituting these 
two values into (70), and performing the integral over r, 
the Fermi golden rule expression without the FES effects 
emerges in the imaginary part of "f x h-hx representing the 
particle loss rate from the QD. 
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To compute the integrand in (70), we represent K a 
in the plane-wave basis over the two-dimensional quasi- 
momentum within the surface band. The circular sym- 
metry of this band reduces the problem significantly as 
the plane wave basis decomposes into a product of angu- 
lar momentum eigenfunctions, e M , and Bessel functions 
J\ (At) for the total momentum k. In this basis, K a can be 
represented as a block-diagonal matrix with each block 
corresponding to an angular momentum eigenvalue I, and 
given by 

h 2 k 2 

hKl a -kk' = o Skk ' 

2m L 

+ j dr rJ,(*r)J,(fc'r) J dz\^ L . s (z)\ 2 V aa (r, z), 

where is the effective mass of the surface band, 
V aa (r, z) is the potential due to the QD in state a (see 
(37)), \ipL;s(z)\ 2 is the planar averaged probability den- 
sity of the surface state with negligible fc-dependence (see 
Fig. 4, and Sec. IV B), Ji(kr) are Bessel functions nor- 
malized to unity over [0, R] such that Ji(kR) = 0, and 
R is a cutoff radius set by discretization of k and is 
much larger than the screening length within the two- 
dimensional surface electron gas. In addition, V aa (r, z) 
is constructed in accordance with the discussion in the 
previous section where the potential outside and inside 
the electrode is given by (67) respectively, and followed 
by Fourier transform back to real space. 

We compute the resulting matrices K l a . kk , over a dis- 
crete set of k, and diagonalize them to compute their 
exponentials in (70), numerically exactly. Multiplying 
by Tf k T^ k , and summing over all k, k! in the discrete set, 
we numerically compute the integrand for a discrete set 
of r. By experimenting with the discretization of k and 
the total number of angular momenta /, we obtained well- 
converged results by using 800 k points over the surface 
band shown in Fig. 4, and setting maximum I to 24. 

The resulting integrand as a function of r generally 
has a slowly decaying tail that prevents a direct applica- 
tion of Fourier transform to obtain 7(w). We follow the 
well-established methods to handle this numerical tech- 
nicality 24,28 ' 65 , and then Fourier transform the resulting 
expression to obtain 7(w). The remaining procedure to 
obtain B^h-xx is straightforward. 

In Fig. 6(a-d) we show b x (t) for two different align- 
ments between the energy level of the lowest electron 
state and the Fermi level of the electrode. Results are 
also shown for two different temperatures, and plots in 
each figure correspond to the absence and presence of 
the Coulomb interaction between the QD and the elec- 
trode. We see from the figure that the tunneling rate 
in the presence of the interaction is always smaller than 
in the absence of this coupling. This effect is the result 
of the FES, which in turn is mainly dominated by the 
AOC function rather than the Mahan exciton contribu- 
tion. We verified this by comparing these results with 
calculations in which the AOC is excluded. Decrease in 



the rate also results from removal of the substrate be- 
cause it eliminates screening of the QD potential. 

We observe in Fig. 6 that the suppression in the tun- 
neling rate increases with temperature. This trend fol- 
lows from the Anderson- Yuval mapping 30 , from which we 
expect the AOC function to exponentially decaying for 
time t > h/kgT and at temperatures much smaller than 
the Fermi temperature of the electron gas. Since the 
present calculation is precisely within this regime, the 
suppression of rate increases with temperature. At tem- 
peratures exceeding the Fermi temperature, the orthogo- 
nality catastrophe would itself become exponentially sup- 
pressed. This regime could be reached with an electrode 
made of lightly doped semiconductor in which the Fermi 
energy can be an order of magnitude smaller. This is also 
the kind of system studied in an experiment by Klcemans 
et. al. . 

In addition, note from Fig. 6(b-c) that the time for 
b x (t) to reach its asymptotic value is much smaller than 
the inverse rate implied by its magnitude. Thus the tun- 
neling process may be modeled accurately as Markovian 
with a constant rate given by b x (t) for t > lps. At tem- 
perature of 10K and below, as shown in Fig. 6(a), non- 
Markov behavior may be expected. The longer time for 
the approach to asymptotic value in this regime is the 
result of sharp increase in the density of vacant states. 

We conclude that for a semi-infinite metallic elec- 
trode, the edge singularity effect is small and quanti- 
tative, which is also in agreement with recent ab-initio 
work on the topic 23 . On the other hand, since the effect 
is proportional to the ratio of the scattering potential to 
the bandwidth of the Fermi lightly doped semicon- 

ductor would be a better system for its observation. 

Let us briefly comment on how this crossover from 
the Markov to non-Markov regime may be accessible for 
observation. Since it occurs within the sub-picosecond 
timescale, the photoluminescence in this regime is com- 
pletely quenched. The crossover is therefore relevant 
mainly to non-linear optical response of the system. A 
possible route to accessing this in nonlinear optics is the 
reduction in bleaching of exciton absorption line due to 
dissociation. This bleaching can be studied as a function 
of delay between pump pulses tuned to the absorption 
frequency. The rise in absorption as a function of the 
delay would then change from an exponential to a non- 
exponential function as the temperature is lowered across 
the crossover, which occurs at approximately 10K in the 
present model. 



E. Non-radiative exciton recombination 

In addition to dissociation, photoluminescence may 
also be quenched by non-radiative recombination (NRR) 
in which the energy is transferred to the electron-hole 
pair excitations in the electrode, rather than being con- 
verted to a photon. In the present geometry, this process 
is also faster than the typical radiative recombination of 
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excitons in InAs, as we now discuss. 

We first quantify the NRR of an exciton and its sensi- 
tivity to various physical properties of the system. With 
reference to the full rate expression in (39), we focus on 
M gg:xx so Dgg and D xx enter , both of which are sim- 
ply unity. Thus our theory describes this process by the 
matrix element C gg - xx (t), which depends on d c „ • E(r), 
where d c „ is the conduction-valence band dipole matrix 
element at the band edge. We set its value to 2.15 A 66 . 
The vector E is the total electric field of the transition 
density of the exciton (see Appendix A) , 

E(r) = -V J dsV{r;s)ip n (s)ip n ,(s), 

where the potential V also includes the electronic surface 
polarization response of the QD. 



1 p+oo 

C gg -xx{t) = K2^2\V xg (q)\ 2 / dujX (q, u)S(u-u) gx , t) , 

n q J-oo 

(71) 

where X(cj) is the susceptibility the computation of 
which we described above in Sec. IV C, and 



V xg (q)=d cv - j d?re- i( * r K{r). 



(72) 
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Figure 6: The escape rate of the electron in the 
presence of a hole with (solid red line) and without 
(dashed-dotted line) Coulomb coupling to the electrode. 
The triple-dotted green line corresponds to the rate in 
the absence of the bulk substrate, but with interaction 
with the 2D surface gas retained. Note the suppression 
of tunneling due to the FES. Only the lowest electron 
level and all the hole levels correlated to it are included. 
The effective junction width is 6.25 A. Plots are shown 

for two different temperatures for the Fermi sea of 
surface states, and AE is the energy difference between 

the Fermi level and the lowest electron level. The 
oscillations at the rightmost end of each plot are due to 
the onset of aliasing as the time increases beyond the 
inverse of the largest energy differences imposed by the 
2D momentum discretization of plane wave surface 
states. 



We compute the integrand over a discrete set of (q, uS) 
over the range shown in Figs. 5, and sum over this range 
to compute C gg - tXX (t) at a discrete set of points t. The 
range in q is sufficient due to the fact that the transi- 
tion density of the exciton has a real space spread of 
approximately 4 nm. The range in u) is justified from 
the exponentially suppressed X(q, lj) above the acoustic 
plasmon line in Fig. 5. 

In Fig. 7, we plot C gmxx (t) for various values of I at 
which the multipole expansion of E is truncated. Each 
I > plot is therefore a correction to the often used 
point dipole model, in which only the I = term of the 
transition densities is employed. We see from the figure 
that the converged solution is almost twice as large as 
I = case, and that convergence occurs only beyond 
I = 3. 

As in the case of charge transfer, we may also model 
NRR as a Markov process with a constant rate equal to 
the asymptotic value of C gg - xx (t). This is verified by plots 



of C a 



.(t) in Fig. 7, which show that after the initial 



appearance of an exciton state, the Fermi sea responds 
at the ultrafast timescale of 10 fs, within which C gg - xx (t) 
settles to a constant value. This is much faster than the 
rate implied by the asymptotic value. 

We remark that while NRR is almost 100 times smaller 
than electron tunneling rate in the figures shown, we ex- 
pect from the exponential suppression of tunneling with 
distance, as opposed to a power law dependence of the 
NRR, that the non-radiative decay of an excited QD 
would cross over from a dissociative to a non-dissociative 
channel as its distance from the electrode is increased. 
In the present model, the tunneling rate is reduced by 
a factor of 10 for every extra l.lA of separation so that 
both regimes may be accessed for systems with only sub- 
nanometer differences in the tunnel junction. Thus the 
efficiency of current extraction in a photovoltaic device 
may be heavily impacted by NRR, and a point dipole 
model with uncorrelated electron-hole density suffices to 
estimate this impact within a factor of 2 as can be seen 
from Fig. 8. The correlation between the electron and the 
hole slightly increases the decay rate by shifting the tran- 
sition density towards the metallic surface (see Fig. 9) 



F. Hole cooling and tunneling 

Let us now turn to hole tunneling as well as cooling via 
Coulomb driven energy transfer to pair excitations. The 
crossover between charge and energy transfer as the dom- 
inant decay channel is also relevant here. To show this 
we calculate the matrix elements B gg -hh(t) for tunnel- 
ing, and Chh-h>h'{t) for Coulomb-driven cooling. There 
is a symmetry restriction in the latter so that transitions 
conserve the angular momentum along the surface nor- 
mal of the electrode. Furthermore, due to our approx- 
imation [K x ,Hl] ~ in Sec. IIIB, the computation of 
B gg -hh{t) is obtained by replacing x by g in B xx . hh (t). 
Similarly, the computation of C'hh-h'h' (t) closely paral- 
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Ms that of C gg - xx (t) described in Sec. IIIC. Therefore, 
we have omitted a discussion of computational details in 
this sub-section. 

We plot the results of our calculations for energy trans- 
fer in the presence of a semi-infinite bulk electrode in 
Fig. 10. To understand the dynamical effects of bulk 
screening, we repeated this calculation by keeping the 
static image potential of the QD, but removed the screen- 
ing of electron-electron interaction within the 2D gas. 
Thus we capture only the dynamical effects of the changes 
in the plasmon branch, and plot the resulting rates in 
Fig. 11 (the underlying loss function for this calculation is 
shown in left panel of Fig. 5). The reduction in the rates 
compared to the screened case can be understood by ob- 
serving that only the points in proximity of q ss 0.5 /nm 
and E = 0.2 eV couple to the QD potential. This re- 
gion, which lies entirely inside the incoherent pair exci- 
tation regime is farther from the plasmon branch in the 
absence than in the presence of the substrate screening. 
The resulting loss function is therefore smaller in the un- 
screened than in screened case. Thus the dynamical effect 
of screening is to enhance the rate. 

However, we emphasize that the largest quantitative 
effect of bulk screening is to suppress this rate via the 
instantaneous image potential within the plane of the 
surface electron gas. In Fig. 12 we plot the rates when 
the image potential is completely removed, which corre- 
sponds to the lack of any substrate. Due to the fact that 
the coupling is proportional to the square of the QD po- 
tential, we see a substantial increase, by approximately 
a factor of 250, in the cooling rate. 

We now consider these results in light of the tunneling 
rates of holes, plotted in Fig. 13. In the presence of the 
substrate, tunneling is much faster for s (I = 0) hole 
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Figure 7: Forster rate calculated for increasing order / 
(values indicated in legend) of the multipole expansion 
of the transition density of the exciton, including the 
electrostatic reaction field of the quantum dot. The 
oscillations result from the acoustic plasmon branch in 
the rcnormalized susceptibility of the surface electron 

gas. 
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Figure 8: The Forster decay rate calculated for a 
product state and the fully correlated state for l z = 0. 
The slight increase in the Forster rate for the correlated 
state is the result of shift in the transition density as 
shown in Fig. 9. 
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Figure 9: Contours of transition density of the fully 
correlated exciton state (left) compared to the charge 
density of the product of lowest energy electron and 
hole wavefunctions (right). The thickness gray contour 
lines represents the contour level, ranging from 1 
(maximum) to (minimum). 



states, but approaches the energy transfer rates for p (I — 
1) states. The much smaller tunneling rate for the p state 
is expected from its spatial profile which contains a node 
in the tunnel junction. In the absence of the bulk, the 
energy and charge transfer rates for all states lie within 
factors of two to four and the two processes can thus 
compete. We expect this to hold true even when we 
account for the effect of substrate removal on tunneling. 

To see why, we note that the 2D gas also generates an 
image potential away from its own plane, which would 
have the same effect on the junction potential as the im- 
age potential of a metallic substrate. The formation of 
this image charge may be treated as instantaneous for 
the purpose of tunneling, because as we have seen in the 
results above, the plasmon oscillations are in the regime 
of 10-20 fs, and the dynamical rates settle to constant 
values beyond this timescale. Thus the required charge 
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Figure 10: Cooling rates of hot holes in which energy is 
transferred to the electrode via the Coulomb coupling 
to its dynamical dielectric function. By symmetry the 

transfer conserves l z , as indicated along each curve 
where the change in the principal quantum number n is 
also indicated. The static screening by bulk substrate is 
included both in the response of the surface electron 
gas, and as an image correction to the field of the 
quantum dot. 



re-arrangement in the 2D gas is still much faster than 
the tunneling rates. Also, as shown in Fig 6, the sub- 
strate removal reduces tunneling rates due to the FES, 
but this is still within a factor of 2. The reduced dimen- 
sionality of the electrode thus only lowers the screening 
of the external field within the plane of the 2D gas, and 
has only a small effect on the charge transfer. We con- 
clude that tuning the substrate allows energy transfer to 
be controlled almost independently of charge transfer. 



V. DISCUSSION 

In the previous sections, we have derived a microscopic 
theory and an effective model of dynamics to obtain a 
self-contained framework for studying surface coupled 
quantum dots. In this section, we discuss how the model 
can be applied to analyze a realistic experimental sce- 
nario, and also to analyze the experimental data. We also 
discuss how various physical processes, neglected here for 
brevity, can be described within our framework without 
any further modifications of the theory. 

As an experiment would generally involve a collec- 
tion of quantum dots spread over a region possibly much 
larger than their size, the charge and energy transfer rates 
must be averaged across the collection, and the spatial 
dependence of the electrode density of states. Thus both 
types of rates must be computed for a distribution of 
tunnel junction widths and changes in density of surface 
states across the collection of QDs. Similarly, while the 
optical wavelength is much larger than the size of the 
QD, the total optical response must also account for Spa- 
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Figure 11: Cooling rates for the same states as in 
Fig. 10 but without the substrate induced screening of 
the electron-electron interaction in the surface electron 
gas. The reduction in the cooling rate is almost by a 
factor of 8, and is mainly due to the shifting of the 
plasmon branch up in energy or farther away from the 
part of loss function that couples to the transition 
densities. Line styles are as in Fig. 10. 
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Figure 12: Cooling rates for the same states as in 
Fig. 10 but after completely excluding the effects of the 
bulk substrate. The almost 200-fold increase in the rate 
is due to the removal of image cancellation due to the 
static substrate response. A low dielectric substrate 

supporting a conducting thin film, such as 
metal-on-oxide, would correspond to this case. Line 
styles are as in Fig. 10. 



tial changes in the phase of the waves across the region 
containing the QDs. This could be taken into account by 
including the factor e tk ' Vj in the optical field, where k is 
the wave-vector of the field and r j points to the center 
of of the j th dot. 

Quantum dot arrays may be designed to exploit the 
vast range of timescale of tunneling and the qualita- 
tive changes in exciton decay pathways, which we have 
demonstrated above, in order to control the transport 
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Figure 13: Rate of tunneling into a hole level of the QD 
with a single hole. Plots are shown for lowest energy 
states (principal quantum number n) at two different 
values of l z . The charge distribution of the l z = 1 hole 
state is minimal near the tunnel junction, and therefore, 
it has a much smaller escape rate than the l z — state. 



of energy and charge across these arrays. Our theory 
can be extended straightforwardly to address these sys- 
tems. The Hilbert space of QD states may be expanded 
to include electronic states of the entire array, with ap- 
propriate expansion in the size of the "system density 
matrix" used in our dynamical model. The Hamiltonian, 
Ho would then involve both the intra-dot energy levels 
and inter-dot couplings, for example in a tight-binding 
form. The coupling to electrodes may then acquire addi- 
tional dependence on QD locations, but they could still 
be calculated using the same methodology as outlined in 
this paper. Population dynamics of the density matrix 
of this dot array, computed in the same manner as for 
a single QD above, would then provide the information 
necessary for energy and charge transport studies. 

In a similar vein, the QD density matrix may also 
be extended to include bi-exciton, tri-excitons, and even 
higher order charge complexes. The fundamental struc- 
ture of our theory is a set of levels for the QD, and their 
dynamical coupling to the electrode. From our semi- 
analytical model of the surface polarization in the sur- 
face coupled QD system, the correlation energy of all 
higher-particle states may be determined by expanding 
them in product of single particle basis (though it could 
be impractical beyond the biexcitons). The calculation 
of hybridization and Coulomb interaction with the elec- 
trode then follows the same path as was used for exci- 



tons. The implications of inter-particle correlations on 
their Coulomb coupling to the electrode surface, and the 
energy transfer rates may be studied within this as was 
done above for exciton decay by NRR. We recall that 
the correlations become important when the monopole 
moment of the charge distribution of a state vanishes, 
and therefore does not set the energy scale for Coulomb 
coupling far above all the multipolc contributions. 

We have neglected vibrational levels of the QD in for- 
mulating the present theory. At the simplest level, the 
electron-phonon coupling within the QD provides an ex- 
tra energy transfer channel for the electronic system. The 
discretization of vibrational levels would generally result 
in a non-Markovian time-dependent rate of energy ex- 
change, which may be accounted for via a model spectral 
density of phonons in the QD. This would result in an 
extra energy transfer matrix defined in the same way as 
C a c:a'c'(t) in (65) above, the spectral function of the sus- 
ceptibility in that formula replaced by the model spec- 
tral density of phonons. However, this treatment would 
still neglect the coherent dynamics of the electronic and 
the lattice system. The most extensive model would be a 
coupled system of electronic, photonic, and vibronic den- 
sity matrices. This would allow a full study of phonons, 
excitons, polaritons, and photolumincsccnce. 



VI. CONCLUSION 

In summary, we have presented a theoretical and com- 
putational framework to model the charge kinetics of 
optically excited quantum dots on surfaces. We have 
started from a microscopic construction of the eigenstates 
of the sub-systems, and the couplings between them, and 
employ them in an effective model of dynamics restricted 
to a small subset of these states. The model is ideally 
suited to explore various questions regarding the effect 
of the surface on charge extraction, exciton lifetime, and 
photoluminescence quenching in these systems. We have 
also developed a the dynamical theory taking into ac- 
count the effects of the FES in the electrode response. 

We illustrated the use of our theory by applying it to an 
In As quantum dot above an Au[lll] surface to which it is 
coupled via the Shockley surface states. Both the charge 
and energy transfer processes are essentially Markovian 
in this system, but complex dynamics may result from 
their interplay. We found that the FES can lower the 
exciton dissociation, but does so only by a factor of 2 for 
a Fermi sea that has a bandwidth much larger than the 
screened Coulomb coupling. 

We also discussed the effects of multipole moments of 
excitons on the rate of energy transfer to the electron-hole 
pair excitations in the surface electron gas. The electron- 
hole correlation and its interplay with the image poten- 
tial of the electrode has significant quantitative effects 
on the oscillator strength for energy loss to the plasma 
excitations of the electron gas. This non-radiative decay 
of excitons was found to be about 100 times slower than 
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the dissociative rate for the material and geometry used 
in the calculations. However, with exponential scaling of 
the dissociative rate with the barrier width and height, 
and for larger QDs exhibiting greater correlation effects, 
competition between the two decay pathways can be ex- 
pected. This would also yield rich dynamics. 

Such a Coulomb driven process can also act as a cool- 
ing mechanism for a charged QD, and we have discussed 
in particular the cooling of holes. The energy transfer 
is very sensitive to the excitation spectrum at the elec- 
trode surface, and in particular the coupling of the QD 
to the acoustic surface plasmons. The size of the QD, 
and the energy level spacing, control the strength of its 
dynamical coupling via this mechanism, and we find that 
the shifting and broadening of the plasmon mode due to 
screening by the bulk yields a large (approximately 200- 
fold) reduction in the cooling rate. By arguing that this 
screening does not affect the formation of image poten- 
tials in the junction, we claim that screening properties 
of the substrate can sensitively tune this system between 
energy transfer and charge transfer modes of operations. 



Our work can be applied to model and analyze experi- 
ments on colloidal quantum dots near semiconductor and 
metallic surfaces, and epitaxial quantum dots in multi- 
ple quantum wells. We have emphasized the vast range 
of timescales that can exist in the dynamics of these sys- 
tems, and point to the various cross-overs in the dom- 
inant decay pathways for exciton states as well as the 
charge kinetics. In a subsequent paper, we will develop 
the theory of linear and non-linear optical response of 
these systems. Thus the excitation process itself will be 
studied dynamically in the presence of charge kinetics 
described here. 

Another interesting application of this model that we 
have briefly discussed is the exploration of spatially de- 
pendent coupling between a quantum dot array and an 
electrode to control transport physics in these systems. 
Finally, extensions of the present theory are necessary to 
describe coupling to the driven vibrational modes of the 
quantum dot, and dynamics of photon degrees of freedom 
for exploring the competition between charge and energy 
exchange with photoluminescence and phonons. These 
extensions will be considered in future publications. 



Appendix A: Matrix elements 
1. Hybridization 

For QD states that are exponentially suppressed at the electrode surface, the hybridization may be defined by 
adding and subtracting Hi, to the total Hamiltonian. The decomposition has the advantage of ease of evaluation in 
terms of potential energy differences outside the electrode. Thus for the single electron or hole QD state we write 

T B v l = J druUr)e~ ik - r [H(r) - H L (r)} <p e {r), (Al) 

tl = J druUr)e- ik r [H(r) - H L (r)} <p h (r), (A2) 
T% = T?? (A3) 
Tll = T h Z (A4) 

Since ip^ represents a single electron orbital in the valence band, such that ip* h is the state of the hole, all of the above 
matrix elements are defined with respect to an electron transfer. Here we have shifted the origin so that z = lies 
at the image plane of the electrode. The integrand vanishes on the electrode side (z < 0), and we simplify the above 
formula within the effective mass approximation with a uniform mass, m inside the QD. Since the kinetic energy 
operator acting on <f> a yields E a , and the function u uk (r) varies on a much faster scale than the rest of the functions 
under the integral, we obtain 

Til = J dr(uUz))e-* k - r U a (r)va(r). (A5) 

Here we have replaced u vk (r) by its average, {u* ku (z)), over the planar [111] unit cell of Au lattice. Outside the 
electrode, {u* kv (z)) decays exponentially, but may have additional dependence on z depending on the pseudopotential 
in the DFT calculation. The potential energy U a in (A5), which depends on the QD state, is defined as, 

(vn \ 
1 [E a -V{r)]e(ren D ). 

Here the pseudopotential Ul for the electrode can be extracted from the DFT calculation, the planar average of 
which in our calculations is shown in Fig. 4. The function Q(r G Qd) is unity inside the QD and zero outside. In 
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the presence of cylindrical symmetry, we may write the above integral in terms of a Bessel transform of order equal 
to the angular momentum quantum number, I, of the QD state, 

T%(*) = (-1)' /2 / drrU a {r,z)va{r,z)Ji{kr), 

and write, 

T2 = J dz{uUz))T%{z). (A6) 

Finally, we introduce matrix elements defining the charge transfer processes that form or dissociate an exciton 
state. Using the coefficients, & e h;x, for the expansion of the exciton state in terms of the electron-hole product states, 
we write the matrix elements for charge transfer involving an exciton state in terms of the above matrix elements 
connecting the ground and single particle states, 

T% = "£$ he ;*T%, (A7) 

e 

h 

T% = T%*, (A9) 
T% = T%*. (A10) 



2. Coulomb interaction 



We write the matrix elements in (20) by separating the classical contribution, 

Vn,vk;n'yk' = ^n,uk\n' ,v'k' ~ &vv'$h,k?W*™ a \ (All) 

In any expectation value, the subtracted term multiplies the sum over occupied states, and therefore produces W^™ 9 ^ 
so long as the trace of the statistical matrix for the electrode is normalized to unity. The matrix elements in (All) 
are given by, 

W n , vk , n >y W = \J drW nn ,(r)e-^ k - k ">- r (A12) 
xu*uk( r ) u »k'(r), 

and 

W { Z 9) = \j ds I drW nn ,(r)Q lmg (r;s), 

where Qi mg { r \ s ) is the charge distribution producing an image potential due to a point charge located at s outside 
the surface. In the above equations, W nn ' is the potential due of the "transitions density" of the two QD states, 
tp n (r)(p n > (r). We determine this potential from our semi-analytical solution 34 for the electrostatic potential, which 
yields the coefficients of multipole moments, Qi m { s ) as a function of s. Following the calculation of Qi m we calculate 
the potential of the transition density as follows. When n and n' refer to the same bulk band, we exploit the fast 
scale of Bloch envelope functions and set their overlap to unity, and obtain 

l,m 

When n and n' belong to different bulk bands, the underlying Bloch envelope functions become orthogonal. In this 
case, we must expand the Coulomb potential over the bulk unit cell of the QD. Let x vary over the unit cell, 

V(r — s — x) = V(r — s) — x ■ VV(r — s). 

When this is averaged over the macroscopic envelope functions ip n (s), we replace W nn i(r) by 

W nn ,{r) + d nn , -VWW(r), 
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where d nn > is the dipole moment between the Bloch states between bands n and n' . This is used in calculating 
energy transfer for exciton states in which case the zeroth order term in the field vanishes. When d nn i also vanishes, 
higher order derivatives of the potential must be coupled to the multiple moments of the microscopic Bloch functions. 
Replacing the overlap of Bloch envelope functions by their planar average, 



O u k,uk'(z) = J dr\\u* vk (r)u vk t(r), 



we obtain, 

(A13) 

= \ I ' dzO vk , vk ,{z) j dr^ k '- k y r m m Ar). 

Appendix B: Charge and energy transfer rates 

In this appendix, we collect the necessary mathematical steps to derive the energy and charge transfer rates as 
defined in (49) and (65). We make use of the equivalence of the following two expressions. Let, 

f(t) = -i f drX(r)e^\ (Bl) 



where X(t) is a causal function, and the factor i is for later convenience. The Fourier transform of X(r), which is 
complex, 

X(oj) = X R {uj)+iXi{u), 
satisfies the Kramers-Kronig relations 29 . Then, we may write f(t) entirely in terms of Xj(u>) as, 

f(t) = \ [ dc^Spec [X(u)] S(u, t), (B2) 



2 . 

where we have defined the spectral function Spec LY(w)] = — 2X/(o;), and 

sin[(w — io )t] 2isin 2 [(w — w )^ 



S(u,t) = — 

U> — LUq U! — LOq 

It is useful to note that the limiting behavior of the function S is, 

lim S(ui, t) = —Tr6(u))t, 

t— ¥ OO 

Ldt 2 

lim S(u. t) = —t + i . 

t->o 2 

1. Charge transfer 

The charge transfer matrix can be described by the functions defined in (50), which are explicitly equal to two of 
the four correlation functions in the definition of B ac . ya i c i(t) in (49). To handle all these correlations, and to derive 
this reduction to two functions explicitly, we define here two additional functions, 

a c -£U*V;*) = -m' - 1") lr c ' b (t")T bc {t')% ca m) , 

r* Ci0a ,(*",f;t) = -m' n (T c ' c {t"m a {t)T aa \t')v\) . 

It follows from the above definitions that 

^ c; w (*V;*) = -K'a;co'(t", t'-,t)]*, (B3) 
A-^ c (t" ) i';*) = -[A» 6 V(t" ) t';t)]* ) (B4) 
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Let us introduce the Fourier transforms of the correlation functions as, 



+00 



(uj,t';t) = / dt"T c , c;aa ,(t",t';t)e 



icot" 



The above definitions of the correlation functions and their relationship, (B3) and (B4) imply that 

Spec[r aa/:c , c (uj,t';t)} = Spec [r^,. a , a (-u},t';t)] , 

Spec [A-,. c ,>, *';*)] = Spec [A c 7 ; c a , a (- W , t'; t) 

We now use the basic relation (B2), and exploit the fact that only the Spec[r(w, t'; t)] and Spec[A(w, t'\ t)] are necessary 
in the final expressions. Thus performing the required permutations of states in the subscript, and setting 

laa'-cd (w, t'; t) = Spec [r aa / ;cc / (oj, t'\ t)} , 
Ka'-cd (w, t'; t) = Spec [A aa / ;cc / (w, t'; t)} , 

we obtain (49) in the main text. We now turn to the evaluation of these spectral functions. 

The functions T and A can be expressed in terms of single particle Green functions of the electrode, within the 
interaction picture discussed in the main body. Here we explicitly develop the expressions for T(uj,t';t), since the 
derivation of A(ui,t';t) follows the same form. These Green functions are physically different depending on whether 
they describe propagation of the system under addition or removal of a particle. Following conventional notation we 
use superscript " >n to describe propagation under addition of an electron to the electrode, and " < " for the removal 
of an electron. Thus letting n c be the electron occupation of the QD in state c, we define 



T / r(f" t) 
1 c'c;aa' \ L t b t b ) 



rS C500 /(t",t / ;t) < = n ° + 1 >< = n * - i. 

r^c,«„/ (t",t';t) n' c = n c - l,n' a = n a + 1, 



(B5) 



where T^, c . aa , (t" , t'; t) describes the tunneling of an electron out of the quantum dot and is given by (51) in the main 
text, while 



iK c ,i' f p -iK c (t'-t) iK a (t"-t) 



-iK a ,t''-PHi 



(B6) 



vv'kk 1 



describes the tunneling from the electrode to the quantum dot. We also define A> in a similar fashion to correspond 
to the ordering of S and c. Note that the superscripts on the hybridization matrix elements, e.g. match the 

pairs of subscripts on the correlation function, and these pairs of indices are constrained by the type of correlation 
function, "lesser" or "greater", in accordance with (B5). 

In the case of a quadratic Hamiltonian, this formula can be further simplified using two mathematical relations 
which can be proven using expressions derived by Hirsch 52 . In particular, our Hl and K operators are all quadratic. 
The Hirsch approach allows us to connect the required full traces to the evaluation of an expression based on those 
operators in a single particle basis. In the following, the left hand side gives the ratio of two thermal averages while 
the right hand side is a specific matrix element in a single particle basis of the enclosed series of operators represented 
by matrices in that basis set: 



AA n x AA m y „ AA2 „t „iA$ — SH r, \ 

e e "Cie c'e e > 



vk N mn (-y,-x) \l+ 



iAm y —iAr, x — iA^ —iAo 



i) N mn (-y, -x)J e 



Amy —iArtX — iAri 

"e e J 



A3 c]e iA ?cie iA " x e' A ™ y e-P H L 



(^iAz giA 2 e iA n x ^iA m y g— {3H L ^ 
vk N nm ( X ,y) [l+ ^nX e iA m y e iA 3e iA 2 _ ^ ^ ^ y) j " X jA n 



x iAm y i A-i 

e e J 



(B7) 



(B8) 



In the second of the two preceding equations, N„ m (x, y) = I — N rim (a;, y), and in the the first equation, 

N m „(y, x) = e lK ™ y e lK - x [l + e~ pHL ] ~* e ~ lK ~ x e~ lK ™ y , 



27 



where the central part is given by (56) in the main text. To apply expressions to (B7 and B8), we substitute. 



e iA n x & iA m y _ e iK c ,t' _ ^.K c ,t^K d (t'-t) 
e iA 2 _ e -iK c (t'-t) e iK a (t"-t) 
e iA 3 =e -iK a ,t"_ 



It is also convenient to define a matrix functions, 

sS,(f',(";() = e- iK *' t 'e iK °' t "e- iK « t "T l ? b '(t)e iK ° t ' -I, 
= e ~ lKc,t ' 7o b ' (t)e tKa,t '' e~ lKat '' e lKct ' - I, 

and a function describing the AOC, 

r cl c;aa>{t",t',t) = ^ e i^f e -iJf„(*'-t) e «f.(t"-t) e -«r.,t" e -flir t y 

To apply these expressions to A^ b ™ a ,, we substitute, 

iA.i iA m y iK c t iK„(t' — t) 
e iA 2 _ e —iK b t' e iK b t" 
e iA 3 =e -iK a ,t" _ 

With these definitions, we obtain the following set of functions to describe electron and hole tunneling between the 
QD and the electrode, 

r>, c ., aa ,(t",t';t) = -i@(t'-t")r c , c ., aa ,(t",t',t) (B9) 
x J2 T^Trt?(vk\-N c , c ,(-t',0) [l + S^; a ,(i'.*";i)N c ' c '(-t',0)] _1 [ e - iK ^' e iK ^"] \v'k'), 

w' kh' 

r< c . aa ,(t",t';t) = iQ(t'- t")r c/c , a , a (t'\t',i) (BIO) 
J2 T ^' T ^' <" fe l N C 'c(*",0) [l + sSi;-(t",i';*)N^(t",0)]^ [e- uc .' t "e ue «' t '] |*/fe') , 

tefe' 

and for A 

A> 5 *L< (*' - *", *';*) = - (t, t", t') (Bll) 

x £ 2S^,(i/fc|N ac (t-t',-t)[l + S^ / (t / ,^ [ e -« c "' e «a""] |fcV), 

A< 6 ;L' (*' - t", t'j t) = -i0(t' - t")r cai6a ' (t, t", (B12) 
x X) ^'T^(^|N ra (i-t',t) [l + S^(f',f";t)N ca (f-i',t)] _1 [ e *.' e *.«-V c ^] |^'fe'). 



2. Energy transfer 

a. First order term 



Since we are only dealing with the Coulomb term in this appendix, we will omit the subscript C used in the main 
text for functions related to this interaction. From the definitions of the superoperators, the first order term takes the 
form shown in (59) in the main text. We now impose the condition that the subscripts of V a b correspond to states with 
the same net charge, and that only the monopole terms are significant in the operators K a . Thus we replace K a and 
K' a by their average, and define the time evolution of the electrode density operator under the potential (V a + V^)/2 
as 
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The corrections to this are small, and quantified at the end of this discussion. The first order term now becomes, 



dt 



Pac:a>c>{t) = -iD ac (t) (t ) 



Clearly, each term under the sum is the expectation value of the potential energy as the electrode state evolves under 
the average potentials of the states that coupled by the off-diagonal Coulomb interaction. Thus by using our definition 
of the Hartree energy matrix in (60), 



d Pac;av(f) - -tfcdUW (^(t) + ^Ei^W f Aqq ,(t) 



dt 



dt 



+i5aa>D-i(t)[A c , c (t) + 



d\ogD ac (t) f f 



dt 



[ A c>c (r) 
Jo 



(B13) 



To obtain the corrections beyond the above, we let K aa = (K a + K a /)/2, and SW aa i — (V aa — V a > a ')/2 J and write 



+sw aa ,)t 



V aa ,e UK U aa ,{t)mJl a ,{t)e 



itK a 



where 



Ua*>(t) =l-i dt'e- ttK 6W aa ,(t')e 
Jo 



l\itK aa 



t ,t' 



e -«'*"' 5W aa , (t')e l(t '-*" )i?aa ' 5W aa , (t")e u " Raa ' + 



From this result, we see that the leading correction to the first order term, as approximated by (B13) is given by a 
potential 5W aa ', which is equal to the difference between two electrostatic potentials of equal charge. 



b. Second order term 



We derive the expression for one of the terms in the energy transfer matrix, while the remaining three terms can 
be computed from it. Thus let, 

F ac . a ' c/ (r;t) = -e i ^' +u ^ t x 

r dT(AV c , c {t-T)AV aa ,(t)m)e 



We define 

xXt'a' fa t) = -i0( T ) (AV c , c (t - r)AV aa , (t)m 
Wa' fa *) = (AVc c (t)AV aa , (t - r)m 

These two functions also satisfy the relation 

(xt'c l aa'( T > t )) = -Xa'a;c'c( T ^) 



Spec 



Xc'c-.aa'i^t) 



Spec 



Xc'c-.aa'i-^t) 



Using the expressions for AV , and implying summation over repeated indices, we obtain, 

- r)c v , kl {t - r)ct p (*) C/iV (t)) 
" (4fc(< - T)c v , k ,(t - r)) (cl P (t)c^ p ,(t)) 
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In most circumstances, we may neglect defects in the electrode surfaces and assume a uniform electron gas. Within 
this approximation, the correlation functions must conserve momentum, and the planar Fourier transform of the QD 
potential also becomes a function only of the difference between the two momenta 



V & , vh;cyv = V c , v , ;cv (k -k') = J V c , c {r)e- i{ - k - k '> r dr. 



In this case the number of momenta in the summation reduces, and we obtain 

xtc-aa'i^i) =^2v c ,, I/ . tCtV ,(-q)X uu/; w>(q,L))Va, fi ; a >,v>(q)- 



Here the function, X vll i ]tlll i{q,uj), is the density density correlation function, generalized to include inter-subband 
transitions, 



A, 



dr e l1 



4k-g(* ~ T K'k{t - T)cl p+q (t) Ct ,, p (t) 

The calculation of density-density correlation can be found in most textbooks on solid state physics 29 . Substitution 
of these calculations into the above formula and setting \ = Spec[x + ] yields formula (65) in the main body. 

I 
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